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Abstract. This is the main paper in a sequence in which we give a complete proof of 
the bounded curvature conjecture. More precisely we show that the time of existence 
of a classical solution to the Einstein-vacuum equations depends only on the L'^-norm 
of the curvature and a lower bound on the volume radius of the corresponding initial 
data set. We note that though the result is not optimal with respect to the standard 
scaling of the Einstein equations, it is nevertheless critical with respect to another, more 
subtle, scaling tied to its causal geometry. Indeed, bounds on the curvature is the 
minimum requirement necessary to obtain lower bounds on the radius of injectivity of 
causal boundaries. We note also that, while the first nontrivial improvements for well 
posedness for quasilinear hyperbolic systems in spacetinie dimensions greater than 1 + 1 
(based on Strichartz estimates) were obtained in [2], [3], [48], [49], [19] and optimized 
in [20], [36], the result we present here is the first in which the full structure of the 
quasilinear hyperbolic system, not just its principal part, plays a crucial role. 

To achieve our goals we recast the Einstein vacuum equations as a quasilinear so(3, 1)- 
valued Yang-Mills theory and introduce a Coulomb type gauge condition in which the 
equations exhibit a specific new type of null structure compatible with the quasilinear, 
covariant nature of the equations. To prove the conjecture we formulate and establish 
bilinear and trilinear estimates on rough backgrounds which allow us to make use of 
that crucial structure. These require a careful construction and control of parametrices 
including error bounds which is carried out in [41]-[44], as well as a proof of sharp 
Strichartz estimates for the wave equation on a rough background which is carried out 
in [45]. It is at this level that the null scaling mentioned above makes our problem 
critical. Indeed, any known notion of a parametrix relies in an essential way on the 
eikonal equation, and our space-time possesses, barely, the minimal regularity needed to 
make sense of its solutions. 

1. Introduction 

This is the main in a sequence of papers in which we give a complete proof of the 
bounded curvature conjecture. According to the conjecture the time of existence of 
a classical solution to the Einstein-vacuum equations depends only on the L^-norm of 
the curvature and a lower bound on the volume radius of the corresponding initial data 
set. At a deep level the curvature conjecture concerns the relationship between the 
curvature tensor and the causal geometry of an Einstein vacuum space-time. Thus, though 
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the result is not optimal with respect to the standard scaling of the Einstein equations, 
it is nevertheless critical with respect to a different scaling, which we call null scaling, 
tied to its causal properties. More precisely. L'^ curvature bounds are strictly necessary to 
obtain lower bounds on the radius of injcctivity of causal boundaries. These lower bounds 
turn out to be crucial for the construction of parametrices and derivation of bilinear 
and trilinear spacetime estimates for solutions to scalar wave equations. We note also 
that, while the first nontrivial improvements for well posedness for quasilinear hyperbolic 
systems in spacetime dimensions greater than 1 + 1 (based on Strichartz estimates) were 
obtained in [2], [3], [48], [49], [19] and optimized in [20], [36], the result we present here 
is the first in which the full structure of the quasihnear hyperbolic system, not just its 
principal part, plays a crucial role. 

1.1. Initial value problem. We consider the Einstein vacuum equations (EVE), 



where RICq^ denotes the Ricci curvature tensor of a four dimensional Lorentzian space 
time {Ai, g). An initial data set for (1.1) consists of a three dimensional 3-surface Eq 
together with a Riemannian metric g and a symmetric 2-tensor k verifying the constraint 
equations. 



where the covariant derivative V is defined with respect to the metric g, Rscai is the 
scalar curvature of g, and trk is the trace of k with respect to the metric g. In this work 
we restrict ourselves to asymptotically fiat initial data sets with one end. For a given 
initial data set the Cauchy problem consists in finding a metric g satisfying (1.1) and an 
embedding of Sq in Ai such that the metric induced by g on Sq coincides with g and the 
2-tensor k is the second fundamental form of the hypersurface Eq C The first local 
existence and uniqueness result for (EVE) was estabhshed by Y.C. Bruhat, see [5], with 
the help of wave coordinates which allowed her to cast the Einstein vacuum equations in 
the form of a system of nonlinear wave equations to which one can apply^ the standard 
theory of nonlinear hyperbolic systems. The optimal, classical^ result states the following. 

Theorem 1.1 (Classical local existence [12] [14]). Let {T,o,g,k) be an initial data set 

for the Einstein vacuum equations (1.1). Assume that Eq can be covered by a locally fi- 
nite system of coordinate charts, related to each other by diffeomorphisms, such that 



The original proof in [5] relied on representation formulas, following an approach pioneered by Sobolev, 
see [37]. 

^Based only on energy estimates and classical Sobolev inequalities. 



RiCa^ = 



(1.1) 



V^kij - VMk = 0, 
Rscai-\k\^ + {trky = 0, 



(1.2) 
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{g, k) G i/fpg(So) X Hi^^{T,q) with s > §• Then there exists a unique"^ (up to an isorn- 
etry) globally hyperbolic development {M.,g), verifying (1.1); for which Eq is a Cauchy 
hypersurface^ . 

1.2. Bounded curvature conjecture. The classical exponents s > 5/2 are clearly 

not optimal. By straightforward scaling considerations one might expect to make sense 
of the initial value problem for s > Sc = 3/2, with Sc the natural scaling exponent for 
based Sobolev norms. Note that for s — Sc — 3/2 a local in time existence result, for 
sufficiently small data, would be equivalent to a global result. More precisely any smooth 
initial data, small in the corresponding critical norm, would be globally smooth. Such 
a well-posedness (WP) result would be thus comparable with the so called e- regularity 
results for nonlinear elliptic and parabolic problems, which play such a fundamental role in 
the global regularity properties of general solutions. For quasilinear hyperbolic problems 
critical WP results have only been established in the case of 1 + 1 dimensional systems, 
or spherically symmetric solutions of higher dimensional problems, in which case the L^- 
Sobolev norms can be replaced by bounded variation (BV) type norms^. A particularly 
important example of this type is the critical BV well-posedness result established by 
Christodoulou for spherically symmetric solutions of the Einstein equations coupled with 
a scalar field, see [7]. The result played a crucial role in his celebrated work on the 
Weak Cosmic Censorship for the same model, see [8]. As well known, unfortunately, the 
BV-norms are completely inadequate in higher dimensions; the only norms which can 
propagate the regularity properties of the data are necessarily based. 

The quest for optimal well-posedness in higher dimensions has been one of the major 
themes in non-linear hyperbolic PDE's in the last twenty years. Major advances have 
been made in the particular case of semi-linear wave equations. In the case of geometric 
wave equations such as Wave Maps and Yang- Mills, which possess a well understood null 
structure, well-posedness holds true for all exponents larger than the corresponding crit- 
ical exponent. For example, in the case of Wave Maps defined from the Minkowski space 
M""*"^ to a complete Riemannian manifold, the critical scaling exponents is Sc = n/2 and 
well-posedness is known to hold all the way down to Sc for all dimensions n > 2. This 
critical well-posedness result, for s = n/2, plays a fundamental role in the recent, large 
data, global results of [46] , [39] , [40] and [28] for 2 -|- 1 dimensional wave maps. 

The role played by critical exponents for quasi-linear equations is much less understood. 
The first well posedness results, on any (higher dimensional) quasilinear hyperbolic sys- 
tem, which go beyond the classical Sobolev exponents, obtained in [2], [3], and [48], [49] 

^The original proof in [12], [14] actually requires one more derivative for the uniqueness. The fact that 

uniqueness holds at the same level of regularity than the existence has been obtained in [33] 
^That is any past directed, in-extendable causal curve in M intersects Sq- 

^Recall that the entire theory of shock waves for 1-1-1 systems of conservation laws is based on BV 
norms, which are critical with respect to the scaling of the equations. Note also that these BV norms 
are not, typically, conserved and that Glimm's famous existence result [13] can be interpreted as a global 
well posedness result for initial data with small BV norms. 
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and [19], do not take into account the specific (null) structure of the equations. Yet the 
presence of such structure was crucial in the derivation of the optimal results mentioned 
above, for geometric semilincar equations. In the case of the Einstein equations it is not 
at all clear what such structure should be, if there is one at all. Indeed, the only specific 
structural condition, known for (EVE), discovered in [30] under the name of the weak null 
condition, is not at all adequate for improved well posedness results, see remark 1.3. It is 
known however, see [29], that without such a structure one cannot have well posedeness 
for exponents^ s < 2. Yet (EVE) are of fundamental importance and as such it is not 
unreasonable to expect that such a structure must exist. 

Even assuming such a structure, a result of well-posedness for the Einstein equations 
at, or near, the critical regularity Sc — 3/2 is not only completely out of reach but may in 
fact be wrong. This is due to the presence of a different scaling connected to the geometry 
of boundaries of causal domains. It is because of this more subtle scaling that we need 
at least L^-bounds for the curvature to derive a lower bound on the radius of injectivity 
of null hypersurfaces and thus control their local regularity properties. This imposes a 
crucial obstacle to well posedness below s = 2. Indeed, as we will show in the next 
subsection, any such result would require, crucially, bilinear and even trillinear estimates 
for solutions to wave equations of the form Dg^ = F. Such estimates, however, depend 
on Fourier integral representations, with a phase function u which solves the eikonal 
equation g°'^daudj3U — 0. Thus the much needed bilinear estimates depend, ultimately, 
on the regularity properties of the level hypersurfaces of the phase u which are, of course, 
null. The catastrophic breakdown of the regularity of these null hypersurfaces, in the 
absence of a lower bound for the injectivity radius, would make these Fourier integral 
representations entirely useless. 

These considerations lead one to conclude that, the following conjecture, proposed in 
[18], is most probably sharp in so far as the minimal number of derivatives in is 
concerned: 

Conjecture[Bounded Curvature Conjecture (BCC)] The Einstein- vacuum equa- 
tions admit local Cauchy developments for initial data sets (Eq, g, k) with locally finite 
curvature and locally finite Li^ norm of the first covariant derivatives of k"^ . 

Remark 1.2. It is important to emphasize here that the conjecture should he primarily 
interpreted as a continuation argument for the Einstein equations; that is the space-time 
constructed by evolution from smooth data can be smoothly continued, together with a time 
foliation, as long as the curvature of the foliation and the first covariant derivatives of 
its second fundamental form remain LP'- hounded on the leaves of the foliation. In fact 
the conjecture implies the hreak-down criterion previously obtained in [26] and improved 
in [31], [51]. According to that criterion a vacuum space-time, endowed with a constant 

^Note that the dimension here is n = 3. 

7As we shall see, from the precise theorem stated below, other weaker conditions, such as a lower 
bound on the volume radius, are needed. 
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mean curvature (CMC) foliation Sj, can be extended, together with the foliation, as long 
as the LjL°°(Et) norm of the deformation tensor of the future unit normal to the foliation 
remains bounded. It is straightforward to see, by standard energy estimates, that this 
condition implies bounds for the LfL'^iT^t) norm of the space-time curvature from which 
one can derive bounds for the induced curvature tensor R and the first derivatives of the 
second fundamental form k. Thus, if we can ensure that the time of existence of a space- 
time foliated by depends only on the norms of R and first covariant derivatives of 
k, we can extend the space-time indefinitely. 

In this paper we provide the framework and the main ideas of the proof of the conjecture. 
We rely on the results of [41], [42], [43], [44], [45], which we use here as a black box. A 
summary of the entire proof is given in [27]. 

1.3. Brief history. The conjecture has its roots in the remarkable developments of the 
last twenty years centered around the issue of optimal well-posedness for semilinear wave 
equations. The case of the Einstein equations turns out to be a lot more complicated due 
to the quasilinear character of the equations. To make the discussion more tangible it is 
worthwhile to recall the form of the Einstein vacuum equations in the wave gauge. As- 
suming given coordinates a;°, verifying H^x" = 0, the metric coefficients ga/s — g{da, d^), 
with respect to these coordinates, verify the system of quasilinear wave equations, 

g^''d^d,g^p^F^p{g,dg) (1.3) 

where Fa/3 are quadratic functions of dg, i.e. the derivatives of g with respect to the 
coordinates x"'. In a first approximation we may compare (1.3) with the semilinear wave 
equation, 

□(/) = F((/., a(/)) (1.4) 

with F quadratic in 90. Using standard energy estimates, one can prove an estimate, 
roughly, of the form: 

110(^)11. < l|0(O)||.exp (c\j\mr)h^dr 

The classical exponent s > 3/2 + 1 arises simply from the Sobolev embedding of H^, 
r > 3/2 into To go beyond the classical exponent, see [34], one has to replace 

Sobolev inequalities with Strichartz estimates of, roughly, the following type, 

(^j\mr)\\l^dT^ ' <C (^1190(0)11^1+. + ^*||n0(r)bi+.^ 

where e > can be chosen arbitrarily small. This leads to a gain of 1/2 derivatives, i.e. 
we can prove well-posedness for equations of type (1.4) for any exponent s > 2. 

The same type of improvement in the case of quasilinear equations requires a highly 
non-trivial extension of such estimates for wave operators with non-smooth coefficients. 
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The first improved regularity results for quasilinear wave equations of the type, 

g^"'{(t))di,d,4) = F{4),d(t)) (1.5) 

with g^'^{(j)) a non-linear perturbation of the Minkowski metric m^"^, are due to [2], [3], 
and [48], [49] and [19]. The best known results for equations of type (1.3) were obtained in 
[20] and [36]. According to them one can lower the Sobolev exponent s > 5/2 in Theorem 
1.1 to s > 2. It turns out, see [29], that these results are sharp in the general class of 
quasilinear wave equations of type (1.3). To do better one needs to take into account 
the special structure of the Einstein equations and rely on a class of estimates which go 
beyond Strichartz, namely the so called bilinear estimates^. 

In the case of semihnear wave equations, such as Wave Maps, Maxwell-Klein-Gordon 
and Yang-Mills, the first results which make use of bilinear estimates go back to [15], [16], 
[17]. In the particular case of the Maxwell- Klein- Gordon and Yang-Mills equation the 
main observation was that, after the choice of a special gauge (Goulomb gauge), the most 
dangerous nonlinear terms exhibit a special, null structure for which one can apply the 
bilinear estimates derived in [15] . With the help of these estimates one was able to derive a 
well posedness result, in the flat Minkowski space ]R^+'^, for the exponent s — Sc + 1/2 = 1, 
where Sc = 1/2 is the critical Sobolev exponent in that casc^. 

To carry out a similar program in the case of the Einstein equations one would need, 
at the very least, the following crucial ingredients: 

A. Provide a coordinate condition, relative to which the Einstein vacuum equations 

verifies an appropriate version of the null condition. 

B. Provide an appropriate geometric framework for deriving bilinear estimates for the 
null quadratic terms appearing in the previous step. 

C. Construct an effective progressive wave representation (parametrix) for solu- 
tions to the scalar linear wave equation = F , derive appropriate hounds for 
both the parametrix and the corresponding error term E = F — Dg^F and use 
them to derive the desired bilinear estimates. 

As it turns out, the proof of several bilinear estimates of Step B reduces to the proof 
of sharp L^{M.) Strichartz estimates for a localized version of the parametrix of step G. 
Thus we will also need the following fourth ingredient. 

D. Prove sharp L^{A4) Strichartz estimates for a localized version of the parametrix 
of step C. 

Note that the last three steps are to be implemented using only hypothetical bounds 
for the space-time curvature tensor, consistent with the conjectured result. To start with. 



Note that no such result, i.e. weh-posedness for s = 2, is presently known for either scalar equations 
of the form (1.5) or systems of the form (1.3). 

^This corresponds precisely to the s = 2 exponent in the case of the Einstein- vacuum equations 
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it is not at all clear what should be the correct coordinate condition, or even if there is 
one for that matter. 

Remark 1.3. As mentioned above, the only known structural condition related to the 
classical null condition, called the weak null condition [30], tied to wave coordinates, fails 
the regularity test. Indeed, the following simple system in Minkowski space verifies the 
weak null condition and yet, according to [29], it is ill posed for s — 2. 

□0 = 0, □t/' = ■ A0. 

Other coordinate conditions, such as spatial harmonic^^, also do not seem to work. 

We rely instead on a Coulomb type condition, for orthonormal frames, adapted to a 
maximal foliation. Such a gauge condition appears naturally if we adopt a Yang-Mills 
description of the Einstein field equations using Cartan's formalism of moving frames^^, 
see [6] . It is important to note that it is not at all a priori clear that such a choice would do 
the job. Indeed, the null form nature of the Yang-Mills equations in the Coulomb gauge 
is only revealed once we commute the resulting equations with the projection operator 
V on the divergence free vectorfields. Such an operation is natural in that case, since V 
commutes with the flat d'Alembertian. In the case of the Einstein equations, however, the 
corresponding commutator term [Dg,?^] generates^^ a whole host of new terms and it is 
quite a miracle that they can all be treated by an extended version of bilinear estimates. 
At an even more fundamental level, the flat Yang-Mills equations possess natural energy 
estimates based on the time symmetry of the Minkowski space. There are no such timelike 
Killing vectorfield in curved space. We have to rely instead on the future unit normal to 
the maximal foliation whose deformation tensor is non-trivial. This leads to another 
class of nonlinear terms which have to be treated by a novel trilinear estimate. 

We will make more comments concerning the implementations of all four ingredients 
later on, in the section 2.4. 

Remark 1.4. In addition to the ingredients mentioned above, we also need a mechanism 
of reducing the proof of the conjecture to small data, in an appropriate sense. Indeed, 
even in the fiat case, the Coulomb gauge condition cannot be globally imposed for large 
data. In fact [17] relied on a cumbersome technical device based on local Coulomb gauges, 
defined on domain of dependence of small balls. Here we rely instead on a variant of the 
gluing construction of [10], [11], see section 2.3. 



Maximal foliation together with spatial harmonic coordinates on the leaves of the foliation would be 
the coordinate condition closest in spirit to the Coulomb gauge. 

^^We would like to thank L. Anderson for pointing out to us the possibility of using such a formalism 
as a potential bridge to [16] . 

^^Note also that additional error terms are generated by projecting the equations on the components 
of the frame. 
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2. Statement of the main results 

2.1. Maximal foliations. In this section, we recall some well-known facts about maxi- 
mal foliations (see for example the introduction in [9]). We assume the space-time {Ai, g) 
to be fohated by the level surfaces of a time function t. Let T denote the unit normal 

to Et, and let k the the second fundamental form of E^, i.e. kab = —g(DaT,eii), where 
Ca, a = 1, 2, 3 denotes an arbitrary frame on E^ and D^T — De^T. We assume that the 
E^ foliation is maximal, i.e. we have: 

tTgk = (2.1) 

where g is the induced metric on Ej. The constraint equations on E( for a maximal 
foliation are given by: 

V"A;„6 = 0, (2.2) 
where V denotes the induced covariant derivative on Ej, and 

Rscai = \k\^. (2.3) 

Also, we denote by n the lapse of the t-foliation, i.e. — — g(Di, Di). n satisfies the 
following elliptic equation on E^: 

An^n\k\\ (2.4) 
Finally, we recall the structure equations of the maximal foliation: 

Wokab = RaObO - n^^VaVftn - kackb (2.5) 



^ahc - ^bkac = R-cOafe (2-6) 

and: 

Rab — kack'' b = ^aObO- (2-7) 

2.2. Main Theorem. We recall below the definition of the volume radius on a general 
Riemannian manifold M. 

Definition 2.1. Let Brip) denote the geodesic ball of center p and radius r. The volume 
radius r^oiip, r) at a point p & M and scales <r is defined by 



rvoi{p,r) = inf 



\Br'{p)\ 



r'<r r"^ 



with \Br\ the volume of relative to the metric on M. The volume radius ryoi{M,r) of 
M on scales < r is the infimum ofr^oiip, r) over all points p e M. 

Our main result is the following: 
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Theorem 2.2 (Main theorem). Let (A^,g) an asymptotically flat solution to the Einstein 
vacuum equations (1.1) together with a maximal foliation by space-like hypersurfaces 
defined as level hypersurfaces of a time function t. Assume that the initial slice (Eg, g, k) 
is such that the Ricci curvature Ric G L^(Eo), V/c G L^(Eo), and Eq has a strictly positive 
volume radius on scales < 1, i.e. r^oii^o, 1) > 0. Then, 

(1) regularity. There exists a time 

T = T{\\Ric\\mEo): ||VA;|U2(s„), r,„KSo, 1)) > 
and a constant 

C = C(||i?zc|U2(s„), ||VA;|U2(2o),r,,KSo, 1)) > 
such that the following control holds on <t <T: 

||R|U-^,L2(s,) < C, ||Vfc|Ujj^jL2(s,) < C and ^M^r^oii^tA) > ^■ 

(2) Higher regularity. Within the same time interval as in part (1) we also have 
the higher derivative estimates^^, 



|a|<m Kl<"^ 



D(«)R|Uc»_^^i2(so < Cm Yl l|V(^)i^^c|U2(2„) + ||V« 



(2.8) 



where Cm depends only on the previous C and m. 



Remark 2.3. Since the core of the main theorem is local in nature we do not need to 
be very precise here with our asymptotic flatness assumption. We may thus assume the 
existence of a coordinate system at infinity, relative to which the metric has two derivatives 
bounded in L"^ , with appropriate asymptotic decay. Note that such bounds could be deduced 
from weighted bounds assumptions for Ric and Vk. 

Remark 2.4. Note that the dependence on ||i?zc||L2(Sy), || VA:||/^2(S|,) in the main theorem 
can he replaced by dependence on ||R||l2(So) where R denotes the space-time curvature 
tensor^^. Indeed this follows from the following well known estimate (see section 8 in 

[26];. 

|VA;|2 + ^|fc|^ < / |R|^ (2.9) 
So 4 

and the Gauss equation relating Ric to R. 



Assuming that the initial has more regularity so that the right-hand side of (2.8) makes sense. 
^^Here and in what follows the notations i?, R will stand for the Riemann curvature tensors of Ef and 
M., while Ric, Ric and RscaU^scai will denote the corresponding Ricci and scalar curvatures. 
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2.3. Reduction to small initial data. We first need an appropriate covering of Sq 
by harmonic coordinates. This is obtained using the following general result based on 
Cheeger-Gromov convergence of Riemannian manifolds. 

Theorem 2.5 ([1] or Theorem 5.4 in [32]). Given Ci > 0, C2 > 0, C3 > 0, there exists Tq > 

such that any 3-dim,ensional, complete, Riemannian manifold {M,g) with \\Ric\\L2(^M) < Ci 
and volume radius at scales < 1 bounded from below by c^, i.e. r„o;(-^)l) ^ ^2, verifies 
the following property: 

Every geodesic ball Br{p) with p & M and r < Tq admits a system of harmonic coordi- 
nates X — {xi^X'i^x^) relative to which we have 

(1 + C3)-'5,,- < g,j < (1 + C3)5,,-, (2.10) 

and 

r I la'^ijf v^rfx < C3. (2.11) 

JbAv) 



We consider e > which will be chosen as a small universal constant. We apply theorem 
2.5 to the Riemannian manifold Sq. Then, there exists a constant: 

ro = ro(||Ric||i2(s„), || V/i:||l2(s„), r^oz(So, l),e) > 

such that every geodesic ball Br{p) with p G Sq and r < ro admits a system of harmonic 
coordinates x = {xi,X2,X3) relative to which wc have: 

(1 + e)-^6ij < gij < (1 + e)5ij, 

and 

I Ic)^ 9i3? V\9\dx < e. 

JBr{p) 

Now, by the asymptotic flatness of Sq, the complement of its end can be covered by the 
union of a finite number of geodesic balls of radius Tq, where the number Nq of geodesic 
balls required only depends on tq. In particular, it is therefore enough to obtain the control 
of R, k and r^oii^t-, 1) of Theorem 2.2 when one restricts to the domain of dependence of 
one such ball. Let us denote this ball by B^q. Next, we rescale the metric of this geodesic 
ball by: 

( \ 

gx{t, = g{^t, Xx), X = mm I J— ^ , p , roe 1 > 0. 

Let^^ Rx, k\ and B^^ be the rescaled versions of R, k and B^q. Then, in view of our choice 
for A, we have: 

ll-RA|U2(B>'g) = V^||-R|U2(B.(,) < e, 
I|VA;a||l2(ba ) = VA||VA;||l2(b^^) < e. 



^^Since in what follows there is no danger to confuse the Ricci curvature Ric with the scalar curvature 
R we use the short hand R to denote the full curvature tensor Ric. 
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and 

Note that B^^ is the rescaled version of 3^^. Thus, it is a geodesic ball for gx of radius 

^ > i > 1. Now, considering gxonQ <t <lis equivalent to considering on < t < A. 
Thus, since Tq, Nq and A depend only on ||-R||l2(Eo)) || VA;||l2(So), r^oii^Oi 1) and e. Theorem 
2.2 is equivalent to the following theorem: 

Theorem 2.6 (Main theorem, version 2). Let (A1,g) an asymptotically flat solution 
to the Einstein vacuum equations (1.1) together with a maximal foliation by space-like 
hypersurfaces defined as level hypersurfaces of a time function t. Let B a geodesic ball 
of radius one in Sq, and let D its domain of dependence. Assume that the initial slice 
(Eoj^j/c) is such that: 

||^||l2(b) < e, ||VA;||l2(b) < e andr^oiiB,!) > ^. 

Let Bf = D n "Ef the slice of D at time t. Then: 

(1) regularity. There exists a small universal constant eo > such that if < 
€ < eo, then the following control holds on <t < 1: 

< e, \\Vk\\L^^^^L2^Bt) ^ e and ^M^r^oi{Bt, 1) > |. 

(2) Higher regularity. The following control holds on < t < 1: 

J2 l|D(")R|Uj^^^jL2(B,) < ||V«i?zc|U2(s) + ||V«V/c|U2(5). (2.12) 

|a|<m 

Notation: In the statement of Theorem 2.6, and in the rest of the paper, the notation 
fi ^ /2 for two real positive scalars /i,/2 means that there exists a universal constant 
C > such that: 

fi < C/2. 



Theorem 2.6 is not yet in suitable form for our proof since some of our constructions 
will be global in space and may not be carried out on a subregion B of Eq. Thus, we 
glue a smooth asymptotically flat solution of the constraint equations (1.2) outside of B, 
where the gluing takes place in an annulus just outside B. This can be achieved using the 
construction in [10], [11]. We finally get an asymptotically flat solution to the constraint 
equations, defined everywhere on Eq, which agrees with our original data set {T,Q,g,k) 
inside B. We still denote this data set by (Eq, g, k). It satisfies the bounds: 

||i?||L2(So) < 2e, ||VA;|U2(So) < 2e and r^oii^o, 1) > \. 



12 



SERGIU KLAINERMAN, IGOR RODNIANSKI, AND JEREMIE SZEFTEL 



Remark 2.7. Notice that the gluing process in [10] -[11] requires the kernel of a certain 
linearized operator to he trivial. This is achieved by conveniently choosing the asymptoti- 
cally flat solution to (1.2) that is glued outside of B to our original data set. This choice 
is always possible since the metrics for which the kernel is nontrivial are non generic (see 

W- 

Remark 2.8. Assuming only L/^ hounds on R and Vfc is not enough to carry out the 
construction in the above mentioned results. However, the problem solved there remains 
subcritical at our desired level of regularity and thus we believe that a closer look at the 
construction in [10] -[11], or an alternative construction, should be able to provide the 
desired result. This is an open problem. 

Remcirk 2.9. Since ||/^|||4(So) ^ ||-R^c||2,2 we deduce that ||/c||l2(_b) ^ e^^^ on the geodesic 
ball B of radius one. Furthermore, asymptotic flatness is compatible with a decay o/|x|~^ 
at infinity, and in particular with k in L^(Eo). So we may assume that the gluing process 
is such that the resulting k satisfies: 

Finally, we have reduced Theorem 2.2 to the case of a small initial data set: 

Theorem 2.10 (Main theorem, version 3). Let (A^,g) an asymptotically flat solution 
to the Einstein vacuum equations (1.1) together with a maximal foliation by space-like 
hypersurfaces defined as level hypersurfaces of a time function t. Assume that the 
initial slice {T:o,g,k) is such that: 

\\R\\l^Eo) < e, II^IU2(So) + ||VA;||i2(So) < e and r^oii^o, 1) > ^■ 

Then: 

(1) regularity. There exists a small universal constant eo > such that if < 
e < eo, the following control holds on <t < 1: 

1 

IIR-IIl-^jL^cs*) < e, ||A;|Uc^_^jL2(s,) + \\Vk\\L^^^^L^i:t) < e and ^M^r^^i^t, 1) > ^■ 

(2) Higher regularity. The following estimates hold on < t < 1 and any m > 0: 

J2 l|D(")R|Uf^^^,L2(s,) < J2 l|V»i^^c|U2(Eo) + ||V»VA;|U2(s„). (2.13) 

|a|<m |'i|<m 

The rest of this paper is devoted to the proof of Theorem 2.10. Note that we will 
concentrate mainly on part (1) of the theorem. The proof of part (2) - which concerns 
the propagation of higher regularity - follows exactly the same steps as the proof of part 
(1) and is sketched in section 13. 
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2.4. Strategy of the proof. The proof of Theorem 2.10 consists in four steps. 

Step A (Yang-Mills formalism) We first cast the Einstein-vacuum equations within 
a Yang-Mills formalism. This relies on the Cartan formalism of moving frames. The 
idea is to give up on a choice of coordinates and express instead the Einstein vacuum 
equations in terms of the connection 1-forms associated to moving orthonormal frames, 
i.e. vectorfields Cq,, which verify, 

g(ea, efs) = niap = diag(-l, 1, 1, 1). 

The connection 1-forms (they are to be interpreted as 1-forms with respect to the external 
index fj, with values in the Lie algebra of so(3, 1)), defined by the formulas, 

(A^)a^ = g(D^e^,eJ (2.14) 

verify the equations, 

D^F^, + [A^F^,] = (2.15) 
where, denoting (F^^)^^^ := Ra0nu, 

(F^,)„^ = (D^A, - D,A^ - [A^, A,])^^. (2.16) 

In other words wc can interpret the curvature tensor as the curvature of the so(3, l)-valued 
connection 1-form A. Note also that the covariant derivatives arc taken only with respect 
to the external indices n^v and do not affect the internal indices a,/3. We can rewrite 
(2.15) in the form, 

□gA, - D,(D^A^) = J,(A, DA) (2.17) 

where, 

J, = D'^([A^,A,])-[A^,F^,]. 

Observe that the equations (2.15)-(2.16) look just hke the Yang-Mills equations on a 
fixed Lorcntzian manifold (A^.g) except, of course, that in our case A and g are not 
independent but connected rather by (2.14), reflecting the quasilinear structure of the 
Einstein equations. Just as in the case of [15], which establishes the well-posedness of 
the Yang-Mills equation in Minkowski space in the energy norm (i.e. s = 1), we rely in 
an essential manner on a Coulomb type gauge condition. More precisely, we take eg to 
be the future unit normal to the St foliation and choose 61,62,63 an orthonormal basis 
to Ej, in such a way that we have, essentially (see precise discussion in section 3.2), 
div A — VMj = 0, where A is the spatial component of A. It turns out that Aq satisfies 
an elliptic equation while each component Ai — g(A, Cj), i — 1,2,3 verifies an equation 
of the form. 



□gA = -di{doAo)+AWjAi + AWiAj+ l.o.t. 



(2.18) 
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with l.o.t. denoting nonlinear terms which can be treated by more elementary techniques 
(including non sharp Strichartz estimates). 

Step B (Bilinear and trilinear estimates) To eliminate 9j(9oAo) in (2.18), we need 
to project (2.18) onto divergence free vectorfields with the help of a non-local operator 
which we denote by V. In the case of the flat Yang-Mills equations, treated in [15], this 
leads to an equation of the form, 

UAi = V{A^djAi) + V{A^diAj) + l.o.t. 

where both terms on the right can be handled by bilinear estimates. In our case we 
encounter however three fundamental differences with the flat situation of [15]. 

• To start with the operator V does not commute with Dg. It turns out, fortunately, 
that the terms generated by commutation can still be estimated by an extended 
class of bilinear estimates which includes contractions with the curvature tensor, 
see section 5.4. 

• All energy estimates used in [15] are based on the standard timehke KiUing vector- 
field dt- In our case the corresponding vectorfield Cq = T [ the future unit normal 
to E^) is not Killing. This leads to another class of trilinear error terms which we 
discuss in sections 8 and 5.4. 

• The main difference with [15] is that we now need bilinear and trilinear estimates 
for solutions of wave equations on background metrics which possess only limited 
regularity. 

This last item is a major problem, both conceptually and technically. On the conceptual 
side we need to rely on a more geometric proof of bilinear estimates based on a plane 
wave representation formula-*^^ for solutions of scalar wave equations, 

□g0 = 0. 

The proof of the bihnear estimates rests on the representation formula^^ 

P POO 

(t)f{t,x)= / e'^'^<*'''^ f{Xuj)\^d\duj (2.19) 

where / represents schematically the initial data^*, and where is a solution of the 
eikonal equation^^, 

Z"^d^''udp''u = Q, (2.20) 

^''We follow the proof of the bihnear estimates outhned in [21] which differs substantiahy from that of 
[15] and is reminiscent of the nuU frame space strategy used by Tataru in his fundamental paper [47]. 

^^(2.19) actually corresponds to the representation formula for a half-wave. The full representation 
formula corresponds to the sum of two half-waves (sec section 10) 

^^Hcrc / is in fact at the level of the Fourier transform of the initial data and the norm ||A/||i2(R3) 
corresponds, roughly, to the norm of the data . 

^^In the flat Minkowski space "u{t, x) =t±x- lo. 
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with appropriate initial conditions on Sq and doo the area element of the standard sphere 
in 

Remark 2.11. Note that (2.19) is a parametrix for a scalar wave equation. The lack 

of a good parametrix for a covariant wave equation forces us to develop a strategy based 
on writing the main equation in components relative to a frame, i. e. instead of dealing 
with the tensorial wave equation (2.17) directly, we consider the system of scalar wave 
equations (2.18). Unlike in the flat case, this scalarization procedure produces several 
terms which are potentially dangerous, and it is fortunate that they can still be controlled 
by the use of an extended^'^ class of bilinear estimates. 

Step C (Control of the parametrix) To prove the bilinear and trilinear estimates of 
Step B, we need in particular to control the parametrix at initial time (i.e. restricted to 
the initial slice Eq) 

cPf{0,x)= / e*^""(°'^)/(Acc;)A2rfA(ia; (2.21) 
7s2 Jo 

and the error term corresponding to (2.19) 

r POO 

Ef{t, x) = □g0/(t, x)=i / e'^ (Dg '^u)f{\uj)X^dXdw (2.22) 

i§2 Jo 

i.e. (pf is an exact solution of Dg^ = only in flat space in which case Dg "^u — 0. This 
requires the following four sub steps 

CI Make an appropriate choice for the equation satisfied by '^u{0,x) on Sq, and 
control the geometry of the foliation o/Eq by the level surfaces of '^u{0,x). 

C2 Prove that the parametrix att = given by (2.21) is bounded in JC(L^{M.^), L^(Eo)) 
using the estimates for ^u{0,x) obtained in CI. 

C3 Control the geometry of the foliation of M. given by the level hypersurfaces of "^u. 

C4 Prove that the error term (2.22) satisfies the estimate ||ii^/||^2(-_^) < C||A/||i2(]R3) 
using the estimates for '^u and Dg '^u proved in C3. 

To achieve Step C3 and Step C4, we need, at the very least, to control Dg '^u in L°°. 
This issue was first addressed in the sequence of papers [22]-[24] where an L°° bound 
for Dg '^u was established, depending only on the norm of the curvature flux along 
null hypersurfaces. The proof required an interplay between both geometric and analytic 
techniques and had all the appearances of being sharp, i.e. we don't expect an L°° bound 
for Dg which requires bounds on less than two derivatives in for the mctric^^. 

To obtain the bound for the Fourier integral operator E defined in (2.22), we need, 
of course, to go beyond uniform estimates for Dg '^u. The classical bounds for Fourier 



^^such as contractions between the Riemann curvature tensor and derivatives of solutions of scalar 
wave equations. 
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integral operators of the form (2.22) are not at all economical in terms of the number of 
integration by parts which are needed. In our case the total number of such integration 
by parts is limited by the regularity properties of the function Dg '^u. To get an bound 
for the parametrix at initial time (2.21) and the error term (2.22) within such restrictive 
regularity properties we need, in particular: 

• In Step CI and Step C3, a precise control of derivatives of '^u and Dg "^u with 
respect to both u as well as with respect to various directional derivatives^^. To 
get optimal control we need, in particular, a very careful construction of the initial 
condition for '^u on Sq and then sharp space-time estimates of Ricci coefficients, 
and their derivatives, associated to the foliation induced by '^u. 

• In Step C2 and Step C4, a careful decompositions of the Fourier integral operators 
(2.21) and (2.22) in both A and cu, similar to the first and second dyadic decom- 
position in harmonic analysis, see [38], as well as a third decomposition, which in 
the case of (2.22) is done with respect to the space-time variables relying on the 
geometric Littlewood-Paley theory developed in [24]. 

Below, we make further comments on Steps C1-C4: 

(1) The choice of u{0,x,u) on Eq in Step CI. Let us note that the typical choice 
u{0,x,u) — X ■ uj in a given coordinate system would not work for us, since we 
don't have enough control on the regularity of a given coordinate system within 
our framework. Instead, we need to find a geometric definition of u{0,x,(x!). A 
natural choice would be 

HgU = on Eo 

which by a simple computation turns out to be the following simple variant of the 
minimal surface equation^^ 

Unfortunately, this choice does not allow us to have enough control of the deriva- 
tives of u in the normal direction to the level surfaces of u. This forces us to look 
for an alternate equation for u: 

= 1 — -; r + k: r on En. 

This equation turns out to be parabolic in the normal direction to the level surfaces 
of u, and allows us to obtain the desired regularity in Step CI. On closer inspection 
it is related with the well known mean curvature flow on Eq. 

(2) How to achieve Step C3. The regularity obtained in Step CI, together with null 
transport equations tied to the eikonal equation, elliptic systems of Hodge type, the 

^^Taking into account the different behavior in tangential and transversal directions with respect to 
the level surfaces of '^u. 

^^In the time symmetric case k = 0, this is exactly the minimal surface equation 
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geometric Littlewood-Paley theory of [24], sharp trace theorems, and an extensive 
use of the structure of the Einstein equations, allows us to propagate the regularity 
on So to the space-time, thus achieving Step C3. 

(3) The regularity with respect to u in Steps CI and C3. The regularity with respect 
to X for u is clearly limited as a consequence of the fact that we only assume 
bounds on R. On the other hand, R is independent of the parameter cu, and one 
might infer that u is smooth with respect to uj. Surprisingly, this is not at all the 
case. Indeed, the regularity in x obtained for u in Steps CI and C3 is better in 
directions tangent to the level hypersurfaces of u. Now, the a; derivatives of the 
tangential directions have non zero normal components. Thus, when differentiat- 
ing the structure equations with respect to ou, tangential derivatives to the level 
surfaces of u are transformed in non tangential derivatives which in turn severely 
limits the regularity in u obtained in Steps CI and C3. 

(4) How to achieve Steps C2 and C4- Let us note that the classical arguments for 
proving bounds for Fourier operators are based either on a TT* argument, or 
a T*T argument, which requires several integration by parts either with respect 
to X for T*T, of with respect to (A,cj) for TT*. Both methods would fail by far 
within the regularity for u obtained in Step CI and Step C3. This forces us to 
design a method which allows to take advantage both of the regularity in x and 
UJ. This is achieved using in particular the following ingredients: 

• geometric integrations by parts taking full advantage of the better regularity 
properties in directions tangent to the level hypersurfaces of u, 

• the standard first and second dyadic decomposition in frequency space, with 
respect to both size and angle (see [38]), an additional decomposition in physi- 
cal space relying on the geometric Littlewood-Paley projections of [24] for Step 
C4, as well as another decomposition involving frequency and angle for Step 



Even with these precautions, at several places in the proof, one encounters log- 
divergences which have to be tackled by ad-hoc techniques, taking full advantage 
of the structure of the Einstein equations. 

Step D (Sharp L^(A^) Strichartz estimates) Recall that the parametrix constructed 
in Step C needs also to be used to prove sharp L^{M) Strichartz estimates. Indeed 
the proof of several bihnear estimates of Step B reduces to the proof of sharp L^{M.) 
Strichartz estimates for the parametrix (2.19) with A localized in a dyadic shell. 
More precisely, let j > 0, and let ip a smooth function on supported in 



Let (pfj the parametrix (2.19) with a additional frequency localization A ~ 2^ 



C2. 




(2.23) 
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We will need the sharp^^ L'^{Ai) Strichartz estimate 

UfjLHM) < 2i||V^(2-^A)/|U2(n,3). (2.24) 

The standard procedure for proving^"' (2.24) is based on a TT* argument which reduces 
it to an L°° estimate for an oscillatory integral with a phase involving '^u. This is then 
achieved by the method of stationary phase which requires quite a few integrations by 
parts. In fact the standard argument would require, at the least^^, that the phase function 
u—'^u verifies, 

dt,,u e L°°, e L°°. (2.25) 

This level of regularity is, unfortunately, incompatible with the regularity properties of 
solutions to our eikonal equation (2.20). In fact, based on the estimates for '^u derived 
in step C3, we are only allowed to assume 

dt,.u e L°°, dt,.d^u e L°°. (2.26) 

We are thus are forced to follow an alternative approach^^ to the stationary phase method 
inspired by [35] and [36] . 

2.5. Structure of the paper. The rest of this paper is devoted to the proof of Theorem 
2.10. Here are the main steps. 

• In section 3, we start by describing the Cartan formalism and introduce compatible 
frames, i.e. frames 60,61,62,63 with 60 the future unit normal to the foliation 
and (61,62,63) an orthonormal basis on S^. We choose 61,62,63 such that the 
spatial components A = (741,712,^3) verify the Coulomb condition VMj = 0. We 
then decompose the equations (3.10)-(3.11) relative to the frame. This leads to 
scalar equations for Aq — g(A, cq) and Ai — g(A, e^) of the form (see Proposition 
3.5), 

A^o = 1-o.t. 

□A = -di{doAo)+A^djAi + A^djAi+ l.o.t. 

where l.o.t. denote nonlinear terms for which the specific structure is irrelevant, i.e. 
no bilinear estimates are needed. The entire proof of the bounded conjecture 
is designed to treat the difficult terms AWjAi and AWjAi. 



Note in particular that the corresponding estimate in the flat case is sharp. 
^^Note that the procedure we describe would prove not only (2.24) but the full range of mixed Strichartz 

estimates. 

^^The regularity (2.25) is necessary to make sense of the change of variables involved in the stationary 
phase method. 

^'''Wc refer to the approach based on the overlap estimates for wave packets derived in [35] and [36] 
in the context of Strichartz estimates respectively for C^'^ and H^^'^ metrics. Note however that our 
approach does not require a wave packet decomposition. 



THE BOUNDED CURVATURE CONJECTURE 



19 



To eliminate c?i(c?ov4o) and exhibit the null structure of the term AWjAi we need to 
project the second equation onto divergence free vectorfields. Unhke the fiat case 
of the Yang-Mills equation (see [17]), the projection does not commute with □ 
and we have to be very careful with the commutator terms which it generates. We 
effectively achieve the desired effects of the projection by introducing the quantity 
B — {—A)~^curl A. The main commutation formulas are discussed in section 6 
and proved in the appendix. 

In section 4, we start by deriving various preliminary estimates on the initial slice 
So, discuss an appropriate version of Uhlenbeck's lemma and show how to control 
Aq, a as well as B = {—A)~^curl A, from our initial assumptions on Eq. 
In section 5, we introduce our bootstrap assumptions and describe the principal 
steps in the proof of version 3 our main theorem, i.e. Theorem 2.10. Note that 
we are by no means economical in our choice of bootstrap assumptions. We have 
decided to give a longer list, than strictly necessary, in the hope that it will make 
the proof more transparent. We make, in particular, a list of bilinear, and even 
trilincar and Strichartz bootstrap assumptions which take advantage of the special 
structure of the Einstein equations. The trilinear bootstrap assumption is needed 
in order to derive the crucial estimates for the curvature tensor. The entire 
proof of Theorem 2.10 is summarized in Propositions 5.7 and 5.8 in which all 
the bootstrap assumptions are improved by estimates which depend only on the 
initial data, as well as Proposition 5.9 in which we prove the propagation of higher 
regularity. 

In section 6, we discuss various elliptic estimates on the slices Ej, derive estimates 
for B from the bootstrap assumptions on A, and we show how to derive estimates 
for A from those of B. 

In section 7, we use the bootstrap assumptions to derive L^-spacetime estimates 
for OB and dOB, estimates which are crucial in order to provide a parametrix 
representation for B and prove the bilinear estimates stated in proposition 5.8. It 
is crucial here that all the commutator terms generated in the process continue to 
have the crucial bilinear structure discussed above and thus can be all estimated 
by our bilinear bootstrap assumptions. 

In section 8, we derive energy estimates for the wave equations Dg^ = F, relying 
again on the bootstrap assumptions, in particular the trilinear ones. 
In section 9, we improve on our basic bootstrap assumption, i.e. all bootstrap 
assumptions except the bilinear, trilinear and Strichartz bootstrap assumptions. 
This corresponds to proving Proposition 5.7. 

In section 10, we show how to construct parametric representation formulas for 
solutions to the scalar wave equation Dg^ = F. The main result of the section. 
Theorem 10.7, depends heavily on Theorem 10.3 whose proof requires, essentially, 
all the constructions and proofs of the papers [41]-[44]. Theorem 10.3 is in fact 
the main black box of this paper. 
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• In sections 11 and 12, we improve on our bilinear, trilinear and Strichartz bootstrap 
assumptions. This corresponds to proving Proposition 5.8. Note that we rely on 
sharp L'^{Jli) Strichartz estimates for a parametrix localized in frequency (see 
Proposition 12.1) which are proved in [45]. 

• Finally, we prove the propagation of higher regularity in section 13. This corre- 
sponds to proving Proposition 5.9. 

The rest of this paper is devoted to the proof of Theorem 2.10. Note that we will 
concentrate mainly on part (1) of the theorem. The proof of part (2), which follows 
exactly the same steps as the proof (with some obvious simplifications) of part (1), is 
sketched in section 13. 

Remark 2.12. To re- emphasize that the special structure of the Einstein equations is 
of fundamental importance in deriving our result we would like to stress that the bilinear 
estimates are needed not only to treat the terms of the form AWjAi and A^djAi mentioned 
above ( which are also present in flat space ) hut also to derive energy estimates for solutions 
to Dg^ = F. Moreover, a trilinear estimate is required to get bounds for R. In addition 
to these, a result such as Theorem 10.3 cannot possibly hold true, for metrics g with our 
limited degree of regularity, unless the Einstein equations are satisfied, i.e. Ric{g) = 0. 
Indeed a crucial element of a construction of a parametrix representation for solutions 
to Dg^ = F, guaranteed in Theorem 10.3, is the control and regularity of a family of 
phase functions with level hypersurfaces which are null with respect to g. As mention a 
few times in this introduction, such controls are intimately tied to the null geometry of a 
space-time, e.g. lower lower bounds for the radius of injectivity of null hypersurfaces, and 
would fail, by a lot, for a general Lorentzian metric g. 

Conclusion. Though this result falls short of the crucial goal of finding a scale invari- 
ant well-posedness criterion in GR, it is clearly optimal in terms of all currently available 
ideas and techniques. Indeed, within our current understanding, a better result would 
require enhanced bilinear estimates, which in turn would rely heavily on parametrices. 
On the other hand, parametrices are based on solutions to the eikonal equation whose 
control requires, at least, bounds for the curvatTirc tensor, as can be seen in many 
instances in our work. Thus, if we are to ultimately find a scale invariant well-posedness 
criterion, it is clear that an entirely new circle of ideas is needed. Such a goal is clearly of 
fundamental importance not just to GR, but also to any physically relevant quasilinear 
hyperbolic system. 
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The works of Bahouri- Chemin [2]-[3] and D.Tataru [49] were the first to go below the 
classical Sobolev exponent s = 5/2, for any quasilinear system in higher dimensions. This 
was, at the time, a major psychological and technical breakthrough which opened the 
way for future developments. Another major breakthrough of the period, with direct 
infiuence on our approach to bilinear estimates in curved spacetimes, is D. Tataru's work 
[47] on critical well posedness for Wave Maps, in which null frame spaces were first in- 
troduced. His joint work with H. Smith [36] which, together with [20], is the first to 
reach optimal well-posedness without bilinear estimates, has also infiuenced our approach 
on parametrices and Strichartz estimates. The authors would also like to acknowledge 
fruitful conversations with L. Anderson, and J. Sterbenz. 

3. Einstein vacuum equations as Yang-Mills gauge theory 

3.1. CEirtan formalism. Consider an Einstein vacuum spacetime (A^,g). We denote 
the covariant differentiation by D. Let Bq, be an orthonormal frame on A4, i.e. 

g(ea, ep) = m«/3 = diag(-l, !,...,!). 

Consistent with the Cartan formalism we define the connection 1 form, 

(A)„^(X) = g(Dxe^,e„) (3.1) 

where X is an arbitrary vectorfield in T[Ai). Observe that, 

(AUiX) = -(A)^„(X) 

i.e. the 1 -form A^dx^ takes values in the Lie algebra of so(3, 1). We separate the internal 
indices a, (5 from the external indices /i according to the following notation. 

(A^)„^ := (A)„^((9^) = g(D^e;3, e„) (3.2) 

Recall that the Riemann curvature tensor is defined by 

R(X, r, [/, V) = g(X, [TiuTiv - -DvTiu - "D^uy]^] ) 

with X, F, U, V arbitrary vectorfields in T{M.). Thus, taking U — d^,V — d^,, coordinate 
vector-fields, 

R(e„, e^, a^, d^) = g(e„, D^D^e^ - D^D^e^). 

We write, 

H^ep = {D^ep,ex)ex = (A^)^y3eA 

and, 

D^D.e^ = D^((A,)%e;,) =9^(A,)%e;, + (A,)%D^e;, 
= 9^(A,)%e;, + (A,)%(A^)%e,. 
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Hence, 

R(e«, Cfi, d^, d^) = df,{A^)ai3 - d^{A^)^p + (A^)„ ^{A^)xp - (A^)„ ^{A^)xp. (3.3) 
Thus defining the Lie bracket, 

([A^, A,])„^ = (A^), ^ (A,)^^ - (A,), ^ (A^)^^ (3.4) 

we obtain: 

or, since d^{Ay) - dy{A^,) = D^A,, - D^A^ 

(F^^)«^ = RarJ^,u = (D^A^ - D^A^ - [A^, AJ)^^. (3.5) 

Therefore we can interpret F as the curvature of the connection A. 

Consider now the covariant derivative of the Riemann curvature tensor, 

— Rd 

— {p.Y,,U - RVg(D./a,/5) - Ri^.g(D./^,/5) 

— ^u^iiv)ah — (Ao-)q, ''(F^,^)^^ — (Ao-)/3 iiy)ci& 

— (D<tF^,,)q,/3 + (Act)q, ''(F^i,)^^ — i^fxy)a ^{A„)l,^ 

= (D,F^. + [A„F^.])^^. 

Hence, 

D,R„^^, = ^^to.F^. D,F^, + [A„ F^,] (3.6) 

where we denote by '-^to the covariant derivative on the corresponding vector bundle. 
More precisely if U = U^j^^.-./Ufc is any /c-tensor on M. with values on the Lie algebra of 
so(3, 1), 

(^)D.U = D.U + [A.,U]. (3.7) 

Remark 3.1. Recall that in {A^)af5, a, (5 are called the internal indices, while fi are called 
the external indices. Now, the internal indices will be irrelevant for most of the paper. 
Thus, from now on, we will drop these internal indices, except for rare instances where 
we will need to distinguish between internal indices of the type ij and internal indices of 
the type Oi. 

The Bianchi identities for Ra/3/^^/ take the form 

(^>D,F^, + (^>D^F,, + (^>D,F,^ = 0. (3.8) 

As it is well known the Einstein vacuum equations Ha/s — imply D^Rq^^ui/ = 0. Thus, 
in view of equation (3.6), 

= (^>D^ F,, = D'^F^. + [A^ F^.] (3.9) 

or, in view of (3.5) and the vanishing of the Ricci curvature of g,. 

□A, - D,(D'^A^) = J, (3.10) 
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where 

J, = D^([A^,A,])-[A^,F^,]. (3.11) 

Using again the vanishing of the Ricci curvature it is easy to check, 

D^J^ = 0. (3.12) 

Finally we recall the general formula of transition between two different orthonormal 
frames Sa and on related by, 

e = O^e 

where rrio,^ = O^^^ ^^755 O is a smooth map from M. to the Lorentz group 0(3, 1). 
In other words, raising and lowering indices with respect to m, 

Oc^aO^^ = 5i (3.13) 
Now, (A^)a^ = g(D^e^,ea). Therefore, 

{KU = OZO^(A^),5 + d,{Ol) (3.14) 

3.2. Compatible frames. Recall that our spacetime is assumed to be foliated by the 
level surfaces St of a time function t, which are maximal, i.e. denoting by k the second 
fundamental form of we have, 

tr^fc = (3.15) 

where g is the induced metric on E^. Let us choose e(o) = T, the future unit normal to 
the foliation, and e(j), i = 1,2,3 an orthonormal frame tangent to S^. We call this a 
frame compatible with our foliation. We consider the connection coefficients (3.2) with 
respect to this frame. Thus, in particular, denoting by Aq, respectively ^j, the temporal 
and spatial components of K^j^ 

(A)oj = {Aj)oi^-hj, = 1,2,3 (3.16) 

(^o)oi = -n-^ViU 1,2,2, (3.17) 

where n denotes the lapse of the i-fohation, i.e. = — g(Df, Di). With this notation 
we note that, 

where V is the induced covariant derivative on Ej and, as before, the notation 'S7i{ki)j or 
V'(Aj)oj, is meant to suggest that the covariant differentiation affects only the external 
index i. Recalling from (2.2) that k verifies the constraint equations, 

V%j = 0, 

we derive, 

V\Ai\^^kr{Ai)mr (3.18) 
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Besides the choice of cq wc arc still free to make a choice for the spatial elements of the 
frame ei, 62, 63. In other words we consider frame transformations which keep Cq fixed, i.e 
transformations of the type, 

ei = 01 Cj 

with O in the orthogonal group 0(3). We now have, according to (3.14), 

{Am)ij = 0^iO]{Am)u + dm{0^i)0]5kl 

or, schematically, 

i„ = OAraO-^ + {d^O)0-^ (3.19) 

formula in which we understand that only the spatial internal indices are involved. Wc 
shall use this freedom later to exhibit a frame ei, 62, 63 such that the corresponding con- 
nection A satisfies the Coulomb gauge condition V\Ai)ij = (see Lemma 4.2). 

3.3. Notations. Wc introduce notations used throughout the paper. From now on, we 
use greek indices to denote general indices on M. which do not refer to the particular 
frame (cq, ei, 62, 63). The letters a, 6, c, d will be used to denote general indices on which 
do not refer to the particular frame (61,62,63). Finally, the letters 1, 171,11 will only 
denote indices relative to the frame (61,62,63). Also, recall that D denotes the covariant 
derivative on Ai, while V denotes the induced covariant derivative on Hf Furthermore, 
d will always refer to the derivative of a scalar quantity relative to one component of the 
frame (60, 61, 62, 63), while d will always refer to the derivative of a scalar quantity relative 
to one component of the the frame (61, 62, 63), so that d — {do, d). For example, dA may 
be any term of the form di{Aj), OqA may be any term of the form dQ{Aj), dAo may be 
any term of the form dj{Ao), and dA = {dA, dAo) = {OoAq, dAo, do A, OA). 

We introduce the curl operator curl defined for any so{3, l)-valued triplet {cui, 002,003) 
of functions on as follows: 

{curlco)i=e^ ^'dj{ooi), (3.20) 

where Gjj; is fully antisymmetric and such that €123= 1. We also introduce the divergence 
operator div defined for any so{3, l)-valued tensor A on as follows: 

divA^V^{Ai)^d^{Ai)+A''. (3.21) 

Remcirk 3.2. Since do and dj are not coordinate derivatives, note that the commutators 
[dj,do\ and [dj,di\ do not vanish. Indeed, we have for any scalar function on A4: 

[di, dj](f) = [ei, ej](f) = (DiCj - 0^-6^)0 

= -((Dj6j, 60) - (Dj-6j, 6o))6o(0) + ((Dj6j, ei) - (Dj6i, 6,))6i(0) 

= -{{Ai)oj - {Aj)oi)do(f) + {{Ai)ij - {Aj)ii)di(l), 
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and: 



[di, do\4> = [ci, eo]0 = (DjCo - Doej)0 

= -((DiCo, Co) - (DoCi, eo))eo(0) + ((DiCo, e,) - (DoCj, e;))e/(0) 
= (Ao)oi9o0 + {{Ai)io - {Ao)H)di(j). 

This can be written schematically as: 



for any scalar function (f) on M.. 

Remark 3.3. The term A? in (3.21) corresponds to a quadratic expression in components 
of A, where the particular indices do not matter. In the rest of the paper, we will adopt 
this schematic notation for lower order terms (e.g. terms of the type A^ and ^) where 
the particular indices do not matter. 

Finally, DAq and DAj will always be understood as □ applied to the scalar functionsAo 
and Ai, while will refer to □ acting on the differential form Ao^. Also, ^Aq will 

always refer to ^{Aq). 

3.4. Main equations for [Aq.,A). In what follows we rewrite equations (3.10)-(3.11) 
with respect to the components Aq and A — (^41,^2, A3). To do this we need the following 
simple lemma. 

Lemma 3.4. For any vectorfield X, we have: 



= Adcf) and [9,, 9o]0 = A^0, 



(3.22) 



X"(nA)e, = n{X ■ A) - ■ DaA - {UX) ■ A. 

Taking X = eo in the lemma and noting that. 



(3.23) 



□ eo 



D"DAeo = -D^(A,)o^e, 



■7 



(A,)o^(A^),'^e^ 



as well as 



,28 



D''(A^) = -Do(Ao) +D^(A,) = - [^qAo + (Ao)o '(A)] + [V^(A)]- (3.24) 




On the other hand 



5o(D'^(A^)) 



dlA,-d,{AoV{A), + do{V\A,)) 



'Recall that trfc = 0. 
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Hence, 

(□A)o - ao(D^(A^)) = nA,-do{AoViA,) + d',Ao + do{A,y{A),-do{V\A,)) 

+ AdA + Ad{Ao) + 

= DAo + d^^Ao - doiV\Ai)) + Ad A + Ad^A^) + A\ 
On the other hand we have, by a straightforward computation, for any scalar 0, 

□0= -ao(ao(/)) + A(/) + n-^Vn- V0, (3.25) 
with A denoting the standard Laplace-Beltrami operator on E^. Therefore, 

(□A)o - 9o(D''(A^)) = AAo-9o(V^(A)) + AaA + AMo + A^ 
Finally, recaUing (3.11), we have, 

Jo = -D'^[A^,Ao] + [A^,F^o]=AaA + A9Ao + A^ 
Hence the cq component of (3.10) takes the form, 

AAo-9o(VMi) = ASA + AMo + A^ (3.26) 
According to (3.18) we have, 

V*(^i)oj — ~ki^[Ai)^j. 

We are thus free to impose the Coulomb like gauge condition, 

W\Ai)jk = 0. (3.27) 

In fact we write both (3.18) and (3.27) in the form, 

V'{A,) = A^. (3.28) 

With this choice of gauge equation (3.26) takes the form, 

AAo = AdA + AOAq + A^. (3.29) 

It remains to derive equations for the scalar components Ai, i = 1,2, 3. First we observe, 
in view of (3.24) and (3.28), 

D^Aa = -DqAo + DM, = -^0^0 + VM, + A^ = -OqAq + A\ (3.30) 

Using lemma 3.4 with X ~ e(i), i = 1, 2, 3 we derive, 

□A = (□A),-2(A^),^D,(A^)-D^(A,),^A^-(A,),^(A,)^^A^ 

or, schematically, ignoring signs or numerical constants in front of the quadratic and cubic 
terms: 

DAi = {DA)i + A^djAi + AodA + AdAo + A^. 
RecaUing (3.12) we have, 

(□A), - di(D'^{A^)) = J,. 
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where Jj is the e(j) component of J. Therefore, 

DAi + di{doAo) = ■djAi+ Ji + AodA + AdAo + A\ 

On the other hand, recalhng the definition of J in (3.11), we easily find, 

Ji = A'-diA + [A^,Fji] + AodA + AdAo + A^. 

Therefore, schematicaUy, 

DAi + di{doAo) = A^ -djAi + A^ -diAj + AodA + AdAo + A^ 

We summarize the results of this subsection in the following proposition. 

Proposition 3.5. Consider an orthonormal frame compatible with a maximal foli- 
ation of the space-time M. with connection coefficients A^ defined by (3.2), their decom- 
position A — (^0,^) relative to the same frame Ca, and Coulomb- like condition on the 
frame, 

divA^ A^. 

In such a frame the Einstein-vacuum equations take the form, 

AAo = AdA + AdAo + A\ (3.31) 

nAi + di{doAo) = AWjAi + A^diAj + AodA + AdAo + A^. (3.32) 

Remark 3.6. It is extremely important to our strategy that we have reduced the covariant 
wave equation (3.10) to the system of scalar equations (3.31) (3.32) (see remark 2.11). 

We also record below the following useful computation. 

Lemma 3.7. We have the following symbolic identity: 

carl {curl {A))j = dj{div A) - A{Aj) + AdA. (3.33) 

Proof. To prove (3.33) we write, using the fact that [di,dj] — Ad in view of (3.22), and 
the definition (3.21) of div : 

curl {curl {A)) j = ejn di{eimn dm{An)) 

dl{dm{An)) 
— {Sjm^ln — (^jn(^lm)dl{dm{An)) 

= m{Ai)) - 

= dj {div A) - A{Aj) + AdA. 
which is (3.33). This concludes the proof of the lemma. □ 
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4. Preliminaries 
4.1. The initial slice. By the assumptions of Theorem 2.10, we have: 

II^IU^CEo) < (4.1) 

miL^^o) + l|V/c|U2(Eo) < e, (4.2) 

and: 

r,„^(Eo,l) > ^. (4.3) 

(4.2), (4.1) and (4.3) together with the estimates in [43] (see section 4.4 in that paper) 
yields: 

11^ - l|U-(So) + ||V^n||i2(So) < e. (4.4) 

Also, we record the following Sobolev embeddings and elliptic estimates on Sq that where 
derived under the assumptions (4.1) and (4.3) in [43] (see section 3.5 in that paper). 

Lemma 4.1 (Calculus inequalities on Sq [43]). Assume that (4.1) and (4.3) hold. We 
have on Eq the following Sobolev embedding for any tensor F: 

II^IU«(Eo) < l|VF|U2(So). (4.5) 

Also, we define the operator (— A)~5 acting on tensors on Eq as: 

"+00 

(i) ' 

where F is the Gamma function, and where U{t)F is defined using the heat flow on Eq.' 

(dr - A)U{t)F = 0, U{0)F = F. 

We have the following Bochner estimates: 

||V(-A)-^|U(^2(s„)) < 1 and ||V2(-A)-i£(^2(s,)) < 1, (4.6) 

where £(L^(Eo)) denotes the set of bounded linear operators on L^(Eo). (4.6) together 
with the Sobolev embedding (4.5) yields: 

||(-A)-^F|U.(,„) < llFll^a^^^^. (4.7) 

4.1.1. The Uhlenbeck type lemma. In order to exhibit a frame ei, 62, 63 such that together 
with eo = T" we obtain a connection A satisfying our Coulomb type gauge on the initial 
slice Eq, we will need the following result in the spirit of the Uhlenbeck lemma^^ [50]. 



2 r+oo 

(-A)-iF=^y T-iu{T)FdT, 



29Note that our smallness assumptions on A make the proof of the Lemma much simpler than the 
original result of Uhlenbeck. 
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Lemma 4.2. Let {M,g) a 3 dimensional Riemannian asymptotically flat manifold. Let 
R denote its curvature tensor and r^oiiM, 1) its volume radius on scales < 1. Let A be a 
connection on M corresponding to an orthonormal frame 61,62,63. Assume the following 
bounds: 

\\A\\l^{m) + W^Ml^m) + \\R\\l^{m) < S and r^oi{M, 1) > ^, (4.8) 

where S > is a small enough constant. Assume also that A and VA belong to l/{M). 
Then, there is another connection A on M satisfying he Coulomb gauge condition V\Ai) — 
0, and such that 

WAWl^m) + ||VA|U2(M) < S (4.9) 

Furthermore, ifW'^A belongs to L'^{M), then V'^A belongs to L'^{M). 

Proof. This is a straightforward adaptation, in a simpler situation, of [50]. Note that in 
the new frame 61, 62, 63, defined by 6j = Ol ej, with O in the orthogonal group 0(3), we 
have, 

A^ = OAmO-^ + {dmO)0-\ 

Our Coulomb gauge condition leads to the elliptic equation for O, 

V"'{{dmO)0-^ + OA^O-^) =0, 0-0^ = I, (4.10) 

which, in view of the smallness assumptions and the boundary condition O — )■ 1 at infinity 
along M, admits the unique solution. We leave the remaining details to the reader. □ 

4.1.2. Control of A, Aq and B = A~^curl{A) on the initial slice. Let us first deduce 
from the Uhlenbeck type Lemma 4.2 the existence of a connection A on Eg satisfying 
the Coulomb gauge condition (3.27). In view of Theorem 2.5, the bound on R in L^(So) 
and on Vyoii^o, 1) assumed in Theorem 2.10 yields the existence of a system of harmonic 
coordinates. Furthermore, let (=1,62,63 an orthonormal frame obtained from 0x^,0x2,9x3 
by a standard orthonormalisation procedure, and let A the corresponding connection. 
Then, the estimates of Theorem 2.5 yield the fact that A and belong to L'^{M). 
Together with the estimates (4.1) on R and (4.3) on r^o;(So, 1), and the Uhlenbeck type 
Lemma 4.2, we obtain the existence of a connection A on Eq satisfying the Coulomb gauge 
condition (3.27). 

Next, using the fact that A satisfies the Coulomb gauge (3.27), and using also the 
estimates (4.1) (4.2) and the estimates of Lemma 4.1 on the initial slice Eq, we may 
estimate A, Aq and B — A~^curl (A). We will make use of the following computation. 

Proposition 4.3. We have the following estimate for A, Aq and B — A~^curl (A) on the 
initial slice Eq.' 
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Proof. We estimate separately the components (Ai)jo, {^i)jb (^o)io and (Ao)ii- We start 
with (A)jo- Recall that (^i)oj = hj- Together with (4.2), we obtain: 

ll('9A)o,|U2(So) < ||VA;|U2(So) + ||A2||l2(So) < e + Plli4(So)- (4.11) 
Also, {Ai)ji = g(Diej,ei) = giViCj^ei). A computation similar to (3.3) yields: 

Thus, we have schematically: 

{curlA)ij = R + A^. 
On the other hand, we have from the Coulomb gauge condition: 

divA = A^. 

Using (3.33), we obtain, writing again schematically: 

(AA)ij- = Vi? + Ad A + A\ (4.12) 
which after multiplication by Aij and integration by parts yields: 

-R||l2(Eo) + II^IIl4(S(,) )||M|U.(Eo) + ||A||i4(So) (4.13) 

< (e+ll^lli4(Eo))II^^IU^(Eo) + ll^lli4(Eo), 

where wc used (4.1) in the last inequality. Now, recall (Aj)oo = 0, which together with 
(4.11) and (4.13) yields: 

l|5^l|L2(Eo)<e+P||i4(So). 

Together with the Sobolev embedding (4.5), this implies: 

l|5^l|L^(Eo)<e. (4.14) 

Next, we estimate Vok. Recall (2.5): 
Also recall Gauss equation (2.7): 

R-aOfcO = Rah — '^kcb- 

Thus, we have: 

Vok^ R-n-^V'^n + A^. (4.15) 
(4.1), (4.4), (4.14), (4.15) and the Sobolev embedding (4.5) imply: 

||VoA;|U2(2„) < e. (4.16) 

Now, {Aj)Qi = kij, and thus: 

{doA)oi = VoA; + AAq, 
which together with (4.16), (4.14) and the Sobolev embedding (4.5) yields: 

||(9oA)odlL2(Eo) < \\^ok\\LHSo) + MLHJ:o)\\M\LH^o) (4.17) 
< e + e||9(Ao)|U.(So). 
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Next, we estimate {doA)ij. In view of (3.3), wc have: 

R(ei, Cj, Co, ei) = {doAi)ij - {diAo)ij + AqA. 

Furthermore, we have: 

Rouj = {d,A,)oi - {djA,)oi + A^^dA + A''. 

Using the symmetry of the curvature tensor Rjjoi — ^Uj, we obtain: 

{doAi)ij = dAo + dA + AA, 

which together with (4.14) and the Sobolev embedding (4.5) yields: 

||(5o^)^,-|U^(Eo) < ||5Ao|U.(s„) + ||M|U.(s„) + ||^|U4(s„)||A|U4(s„) (4.18) 
< e+WdAoh^,:,). 

Since Aqq = 0, (4.17) and (4.18) yield: 

||5oA||l^(Eo) < e+\\dAo\\mj:,). (4.19) 

Next, we estimate 9(^40). Recall (3.31): 

AAo = Ad A + AdAo + A\ 

After multiplication by Aq and integration by parts, and together with (4.14), (4.19) and 
and the Sobolev embedding (4.5), this yields: 

ll^^oll L2(So) ~ (I|A||l4(So)||9A||l2(So) + ||A||l4(Eo)||<9Ao|U2(So) + l|A||i4(So))||^o||L4(Eo) 

which implies: 

\\dAo\\m^,)<e. (4.20) 

Together with (4.19), we obtain: 

HAU^j:,) < e. (4.21) 

Finally, we estimate B on the initial slice Eq using the estimates for A (4.14), (4.20) 
and (4.21). This will be done on in Proposition 6.4. Arguing as in Proposition 6.4 for 
t = together with (4.14), (4.20), (4.21), the Sobolev embeddings (4.5) and (4.7) on Eq, 
the Bochner inequality on Eq (4.6), we immediately obtain: 

\\ddB\\mEo) ^ e- 

This concludes the proof of the proposition. □ 
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5. Strategy of the proof of theorem 2.10 

5.1. Classical local existence. We will need the following well-posedness result for the 
Cauchy problem for the Einstein equations (1.1) in the maximal foliation. 

Theorem 5.1 (Well-posedness for the Einstein equation in the maximal foliation). Let 

{T^Q,g,k) be asymptotically flat and satisfying the constraint equations (1.2), with Ric, 
VRic, k, Vk andV'^k in L^(So), and r^oiiJ^Oi 1) > 0. Then, there exists a unique asymp- 
totically flat solution (A^,g) to the Einstein vacuum equations (1.1) corresponding to this 
initial data set. together with a maximal foliation by space-like hypersurfaces defined 
as level hypersurfaces of a time function t. Furthermore, there exists a time 

T, = T,(|| V(')ffic|U2(s„), < / < 1, II V(^U|U2(s„), < J < 2, r,,KSo, 1)) > 

such that the maximal foliation exists for on < t < T^. with a corresponding control in 
L^,j,^]L2(Et) for Ric, VRic, k, Vk and V^k. 

Theorem 5.1 requires two more derivatives both for R and k with respect to the main 
Theorem 2.2. Its proof is standard and relies solely on energy estimates (as opposed to 
Strichartz estimates of bihnear estimates). We refer the reader to [9] chapter 10 for a 
related statement. 

Remcirk 5.2. In the proof of our main theorem, the result above will be used only in the 
context of an extension and continuity arguments (see Step 1 and Step 3 in section 5.4)- 

5.2. Weeikly reguleir null hypersurfaces. We shall be working with null hyper surfaces 
in Ai verifying a set of assumptions, described below. These assumptions will be easily 
verified by the level hyper surfaces Tiu solutions u of the eikonal equation g^'^d^u d^u = 
discussed in section (10). The regularity of the eikonal equation is studied in detail in 
[43]. 

Definition 5.3 (Weakly regular null hypersurfaces). Let % be a null hypersurface with 
future null normal L verifying g(L, T) — —1. Let also N — L — T . We denote by ^ the 
induced connection along the 2-surfaces T-LnT^f We say that T-L is weakly regular provided 
that, 

||DL|U3(^) + ||D7V|U3(^) < 1, (5.1) 
and the following Sobolev embedding holds for any scalar function f on T-L: 

Wfhe^n) < mh^H) + Wmh^n) + WfhHH). (5.2) 

5.3. Main bootstrap assumptions. Let M > 1 a large enough constant to be chosen 
later in terms only of universal constants. By choosing e > sufficiently small, we 
can also ensure Me is small enough. From now on, we assume the following bootstrap 
assumptions hold true on a fixed interval [0,T*], for some < T* < 1. Note that T-L 
denotes an arbitrary weakly regular null hypersurface, with future directed normal L, 
normalized by the condition g{L,T) — —1. 
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• Bootstrap curvature assumptions 

||R||LrL2(E,) < Me. (5.3) 

Also, 

||R-L|U2(^) <Me, (5.4) 

where R • L denotes any component of R such that at least one index is contracted 
with L. 

• Bootstrap assumptions for the connection A. We also assume that there exist 
A = {Ao,A) verifying our Coulomb type condition on [0, T*] , such that, 

||A||LrL^(EO + II^^IUrL^CEO < Me, (5.5) 

and: 

||^o|Uf'L2(St) + ||^A||LfL2(St) + ||A||L2L°°(St) + || |U?°i:3(St) 

+ ll^^^o||,^.,3^,^^ < Me. (5.6) 

Remark 5.4. Together with the estimates in [43] (see section 4-4 that paper), the 
bootstrap assumption (5.3) yields: 

||^||L-L2(St) + ||V/i;||l-l2(so < Me. (5.7) 

Furthermore, the bootstrap assumption (5.4) together with the estimates in [43] (see section 
4-2 in that paper) yields: 

infr,„,(Ei,l) > ^. (5.8) 

In addition we make the following bilinear estimates assumptions for A and R.: 

• Bilinear assumptions I. Assume, 

\\A^djA\\r.2^M) < Mh\ (5.9) 
Also, for B — (— A)~^citrZ {A) (see (5.37) and the accompanying explanations): 

||/i^a,(as)IU.(^)<MV, (5.10) 

and: 

\\R..,^&B\\L2(M)<Mh\ (5.11) 

Finally, for any weakly regular null hypersurface % and any smooth scalar function 
on M, 

\\kj.d^(t>\WiM) < M\sup\\f<p\\LHH), (5.12) 

n 

and 

\\A^dj4>\\L2^M) < M2esup||y0|U2(^), (5.13) 

H 

where the supremum is taken over all null hypersurfaces T-L. 
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• Bilinear assumptions II. We assume, 



\\{-ArkQv{AAmLHM)<M\^ 



(5.14) 



where the bihnear form Qij is given by Qij{4>, ip) = di(j)djip — djcpdiip . Furthermore, 
we also have: 



sharp. Nevertheless, this estimate will be sufficient for the proof as it will only be 
used to deal with lower order terms. 

Finally we also need a trilinear bootstrap assumption. For this we need to introduce the 
Bel-Robinson tensor. 



We conclude this section by showing that the bootstrap assumptions are verified for some 
positive T*. 

Proposition 5.6. The above bootstrap assumptions are verified on < t < T* for a 
sufficiently small T* > 0. 

Proof. The only challenge here is to prove the existence of the desired connection A, all 
other estimates follow trivially from our initial bounds and the local existence theorem 

above, for sufficiently small T*. More precisely we need to exhibit a frame 61,62,63 such 
that, together with 60 = T, we obtain a connection A satisfying our Coulomb type gauge 
on the slice S^. To achieve this we start on Sq with the orthonormal frame 61,62,63, 
discussed in section 4.1^° and transport it to an orthonormal frame on E^, < i < T*, 
according to the equation. 



■'''such that the corresponding connection A verify the Coulomb gauge condition (3.27) and the esti- 
mates of proposition 4.3 




• Trilinear bootstrap assumption. We assume the following. 



(5.18) 




(5.19) 



T>T{ej)^0, ej{0)^ej, 1,2,3. 
Differentiating, we obtain schematically the following transport equation for A: 

T)t{A) = R, i(0) = A. 
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We can then rely on the estimates of the local existence theorem, for sufficiently small T*, 
to derive Lp'^*jL^(Et) bounds for A, dA and d^A. Since all the bounds for A and R are 
controlled from the initial data, for small T* (thus proportional to e), we are in a position 
to apply Uhlenbeck's lemma 4.2 on to produce the desired connection A. Furthermore, 
differentiating (4.10) twice with respect to Dy, and applying standard elliptic estimates, 
we finally obtain the fact that A, dA and d^A are also controlled in Lj^y^jL^(St) in 
conformity with our bootstrap assumptions. □ 

5.4. Proof of the bounded curvature conjecture. In the following two proposi- 
tions, wc state the improvement of our bootstrap assumptions. 

Proposition 5.7. Let us assume that all bootstrap assumptions of the previous section 
hold for < t < T* . If e > is sufficiently small, then the following improved estimates 
hold true onO<t<T*: 

IML^LH^t) ^ e + M^ei + M^e^ (5.20) 
||R • L||i2(^) <e + M2et+MV, (5.21) 
||^||l-l2(s*) + ||aA|U-L2(s,) < e + Mh^ + Mh\ (5.22) 

||^o||L(°°Z,2(St) + ||^^o||Lt°°i2(St) + ll^o||L2L°°(Et) 

+liaAo|Uoci3(s,) + ||aa^||^^^^3^^^^ < c + mM + mV, (5.23) 

Proposition 5.8. Let us assume that all bootstrap assumptions of the previous section 
hold for < t < T* . If e > is sufficiently small, then the following improved estimates 
hold true onO<t<T*: 

WAWjAU.^M) < (5.24) 
\\A^d,{dB)\\r^.^M)<Mh\ (5.25) 

and 

\\R..jod^B\\L^^M)<M''e\ (5.26) 
Also, for any scalar function (p on M., we have: 



and 



\\k^.d^ct>\\L^^M) < Me (^sup||y0|U.(«) + ||a0|Uc«^.(s^) j , (5.27) 
\\A^dj<P\\LHM) < Me (^sup||y0|U2(^) + Wd^L^LHi:,)^ , (5.28) 
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where the supremum is taken over all (weakly regular) null hypersurfaces %. Finally, we 
have: 

||(-A)-^(g,,(A^))||L^(A4) < Mh\ (5.29) 

\\{-/^)-Hd{A')di{A))\\^.^M) < MV. (5.30) 

Also, 



and 



Il?lt(e.) < Me (5.31) 
l|55||L?Lr(so < Me. (5.32) 



M 



< Mh\ (5.33) 



The proof of Proposition 5.7 is postponed to section 9, while the proof of Proposition 
5.8 is postponed to sections 11 and 12. We also need a proposition on the propagation of 
higher regularity. 

Proposition 5.9. Let us assume that the estimates corresponding to all bootstrap as- 
sumptions of the previous section hold for <t <T* with a universal constant M. Then 
for any t G [0,T*) and for e > sufficiently small, the following propagation of higher 
regularity holds: 

||DR|Uooi2(s,) < 2 (||i?zc|U2(s„) + ||Vi?zc|U2(So) + \\k\\L^(i:,) + ||VA;|U2(s„) + || V2A;|U2(s„)) . 

The proof of Proposition 5.9 is postponed to section 13. Next, let us show how Propo- 
sitions 5.6, 5.7, 5.8 and 5.9 imply our main theorem 2.10. We proceed, by the standard 
bootstrap method , along the following steps: 

Step 1. We show that all bootstrap assumptions are verified for a sufficiently small final 
value T*. 

Step2. Assuming that all bootstrap assumptions hold for fixed values of < T* < 1 and 
M sufficiently large we show that, for e > sufficiently small, we may improve on 
the constant M in our bootstrap assumptions. 
Step 3. Using the estimates derived in step 2 we can extend the time of existence T* to 
T* -\- 5 such that all the bootstrap assumptions remain true. 

Now, Step 1 follows from Proposition 5.6. Step 2 follows from Proposition 5.7 and 

Proposition 5.8. In view of Step 2, the estimates corresponding to all bootstrap assump- 
tions of the previous section hold for < t < T* with a universal constant M. Thus 
the conclusion of Proposition 5.9 holds, and arguing as in the proof of Proposition 5.6, 
we obtain Step 3. Thus, the bootstrap assumptions hold on < i < 1 for a universal 
constant M. In particular, this yields together with (5.7): 

||R|U-L2(E0 < e and ||A:|U-l2(eo < e for all < t < 1. (5.34) 
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In view of (5.8), we also obtain the following control on the volume radius: 

^inf^r,,KSt,l)>^. (5.35) 
Furthermore, Proposition 5.9 yields the following propagation of higher regularity 

J2 l|D(")R|Ufj_^,L^(E,) < Cm ||V«Ric|U2(s,) + ||V«VA;|U2(s,) (5.36) 

|a|<m 

where Cm only depends on m. 

Renicirk 5.10. Note that Proposition 5.9 only yields the case m = 1 in (5.36). The fact 
that (5.36) also holds for higher derivatives m>2 follows from the standard propagation 
of regularity for the classical local existence result of Theorem 5.1 and the hound (5.36) 
with m = 1 coming from Proposition 5. 9. 

Finally, (5.34), the control on the volume radius (5.35) and the propagation of higher 
regularity (5.36) yield the conclusion of Theorem 2.10. Together with the reduction to 
small initial data performed in section 2.3, this concludes the proof of the main Theorem 
2.2. 

The rest of the paper deals with the proofs of propositions 5.7, 5.8 and 5.9. The core 
of the proof is to control A, the spatial part of the connection A. As explained in the 
introduction we need to project our equation for the spatial components A onto divergence 
free vectorficlds. This is needed for two reasons, to eliminate the term 9j(9oAo) on the left 
hand side of (3.32) and to obtain, on the right hand side, terms which exhibit the crucial 
null structure we need to implement our proof. Rather than work with the projection V, 
which is too complicated, we rely instead on the new variable, 

B = {-A)-^curl{A) (5.37) 

for which we derive a wave equation. Since we have (see Lemma 6.5): 

A = curl (B) + l.o.t 

it suffices to obtain estimates for B which lead us to an improvement of the bootstrap 
assumption (5.5) on A. In section 7, we derive space-time estimates for HB and its 
derivatives. Proposition 5.7, which does not require a parametrix representation, is proved 
in 9. Proposition 5.8 is proved in sections 11 and 12 based on the representation formula 
of theorem 10.7 derived in section 10. Finally, Proposition 5.9 is proved in section 13. 



6. Simple consequences of the bootstrap assumptions 



In this section, we discuss elliptic estimates on Ej, we derive estimates for B from the 
bootstrap assumptions on A, and we show how to recover A from B. 
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6.1. Sobolev embeddings and elliptic estimates on S^. First, we derive estimates 
for the lapse n on E^. The bootstrap assumption on R (5.3) and the estimate for k (5.7) 
together with the estimates in [43] (see section 4.4 in that paper) yield: 

Wn - + ||Vn||Loo(x) + \\'V^n\\LooL2(j^^) + || V^n||i,ooi3(S4) (6.1) 

+||V(9on)|U.,3(,,) + ||V^n||^^^^3^^^^ + ||V^(9on)||^^^^3^^^ < Me. 

Remark 6.1. Recall from (3.17) that: 

Thus, the estimates (6.1) forn could in principle he deduced from the bootstrap assump- 
tions (5.6) for Aq. However, notice that Vn e L°°(A4) in view of (6.1), while Aq is only 
in LfL°°{'Et) according to (5.6). This improvement for the components (Ao)oi of Aq will 
turn out to be crucial and subtle^^ (see remark 7.5). 

Next, we record the following Sobolev embeddings and elliptic estimates on that 
where derived under the assumptions (5.4) and (5.3) in [43] (see sections 3.5 and 4.2 in 
that paper). 

Lemma 6.2 (Calculus inequalities on [43]). Assume that the assumptions (5.4) and 
(5.3) hold, and assume that the volume radius at scales < 1 on Eq is bounded from below 
by a universal constant. Let 5 > 0. Then, there exists ro{d) > and a finite covering of 
by geodesic balls of radius ro{5) such that each geodesic ball in the covering admits a 
system of harmonic coordinates x — {xi,X2,X3) relative to which we have 

{l + 5y^5ij<gij<{l + 5)5ij, (6.2) 

and 

ro{S) [ \^^g^J\'^Mdx<S. (6.3) 

Furthermore, we have on Tit the following estimates for any tensor F: 

\\F\\lh^.) < l|VF||,3^^^^, (6.4) 

\\F\\LH^t) < l|VF|U2(so, (6.5) 

||-F||L-(St) < II VF||Lp(St) + ||-F||Lf (Et) Vp > 3, (6.6) 

and: 

^^Using the lapse equation An = n|A:p and /c,V/e S L^L^(St), see (5.7), together with the Sobolev 
embedding (6.5) we only deduce k e L^L^CEt) from which An e L^L^CSt)- This would yield V^n e 
L'^L^{Yit), and thus Vn misses to be in L°°{M.) by a log divergence. However, one can overcome this 
loss by exploiting the Besov improvement with respect to the Sobolev embedding (6.5). We refer the 
reader to section 4.4 in [43] for the details. 
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Finally, we define the operator (— A)~^ acting on tensors on as: 

1 r + oo 

(-A)-^F=-^ / r-tc/(r)Fdr, 

where T is the Gamma function, and where U{t)F is defined using the heat flow on S^; 

{dr - A)U{t)F = 0, U{0)F = F. 
We have the following Bochner estimates: 

||V(-A)-5||^(^.(s,)) < 1 and ||V^(-A)-i£(^2(s,)) < 1, (6.8) 

where denotes the set of bounded linear operators on L^(Et). (6.8) together 

with the Sobolev embedding (6.5) yields: 

||(-A)-iF|U.(^,)<||F||^a^^^^. (6.9) 

Remark 6.3. Note that d^f = V^/ + Adf for any scalar function f on E^. Thus, in 
view of the bootstrap assumption (5.5) for A, we may replace with in the Bochner 
inequality (6.8) when applied to a scalar function. 

6.2. Elliptic estimates for B. Here we derive estimates for B using tlie bootstrap 
assumptions (5.5) (5.6) for A and Aq. 

Proposition 6.4. Let Bi = {—A)~^{curl{A)i). Then, we have: 

l|5(50llLrL^(Eo + l|5'(i^^)IUrL^(s.) + \\didoiB.MLrLH^.) < Me. 

Proof. Using tlie Boclmer inequality on (6.8) together with Remark 6.3, and from the 

bootstrap assumption (5.5) on A, we have: 

\\d{B,)\\L^LH^^) + 11^2(5,) |U,="L^(EO < WAhrLH^^) + \\dA\\L^LH^^) < Me. (6.10) 

Next, we estimate d{do{Bi)). In view of the definition of B, we have: 

do{Bi) = {-A)~\curl (doiA))) + [do, i-A)-']curl (A) + i-A)-\[do, curl]A) 

= {-A)-\curl{do{A))) {~A)-^[do,A]{-A)-\url{A) + {-A)-\[do,curl]A) 

= i-Ar'icurlidoiA))) - i-Ar'[do,A]B + {-A)-\[do,curl]A). 

Thus, in view of the bootstrap assumption (5.5) for A, the Bochner inequality on (6.8) 
and the Sobolev embedding on E^ (6.9), we have: 

||59o(5.)|U^^L2(s,) (6.11) 

< \\d{-A)-\curl{do{A)))h^LHj:,) + \\d{-A)-'[do, A]BU^r.Hj:,) 
+ \\d{-A)-\[do,curl]A)\\L^LH^^) 

< \\doA\\LooL^^,) + \\{-A)-^do,A]BU^LHi:,) + \^^^ 

< Me+ ||(-A)-i[9o, A]S|Uc«^2(s,) + ||(-A)-|([9o, ci.W]A)|Uc«^2(s,) 

< Me+\\[do,A]B\\^^j^^^^ + \\[do,curl]A\\^^j^^^^. 
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Next, we estimate the right-hand side of (6.11). Recall the commutator formula (C.4) 

[^0, A]{Bi) = -2k^'VaVi,iBi) + 2n-'VbnVi,idoiBi)) + n-'AndoiBi) - 2n-^W ank'^'W b{Bi) . 

Together with the Sobolev embedding on (6.5), the bootstrap assumption (5.5) for A, 
and the estimate (6.10) for Bi, this yields: 

ll[^o,A](5.)ll,.,|(,, (6.12) 

^ ll^l|L-L3(s,)||<92(5i)||ioo^2(s,) + ||Vn||iooj^3(s,)||a((9o(5,)) 11^00^2(2,) 

+ l|A^lliooii(S,)ll^0(A)||L-L6(S,) + ||Vn||Loci2(s,)||A;||iooi6(S,)||a(Bi)||ioci6(S,) 

< MV + Me||a(9o(S,))|Uoo^2(s,). 

Next, we estimate the last term in the right-hand side of (6.11). In view of the commuta- 
tor formulas (C.3) and (3.22), and in view of the definition of curl , we have schematically: 

[^0, curl ]A = kVA + n-^VndoA + Ad A = A^A, 

which together with the bootstrap assumption (5.5) for A yields: 

\\[do.curl]A\\^^^.^^^^ < ||A|Uoc^3(s,)||aA|Uc«^2(s^) < M\\ (6.13) 

Finally, (6.11)-(6.13) imply: 

WddoWU^L^i^,) < Me + \\[do,A]B\\^^j^^^^ + \\[do,curl]A\\^^^i^^^^ 

< Me + Me\\d{do{B,mL^LHi:,) 

which yields: 

||5ao(S0llwEO <Me. 
Together with (6.10), this concludes the proof of the proposition. □ 

6.3. A decomposition for A. Recall that B = {—A)~^{curl (A)). We show how to 
recover A from B: 

Lemma 6.5. We have the following estimate: 

A = curl (B) + E 

where E satisfies: 

Proof. In view of Lemma 3.7, we have: 

A = {-A)-^curl {curl (A) + {-A)-\AdA + A^). 

This yields: 

A = curl{-A)-\curl{A)) + [{-A)-\curl]curl{A) + {-A)-\AdA + A^) 
= curl (B) - (-A)-^[A, curl]{-A)-'curl (A) + {-A)-\AdA + A^) 
= curl{B) - {-A)-^[A,curl]B + {-A)-\AdA + A^), 
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which imphes: 

E = -{-A)~'[A,curl]B + {-A)-\AdA + A^). 

Now, we have 

[A, a]0 = Rdcp + dAd(l) + A 

for any scalar function (j) in where the curvature tensor R on is related to R through 
the Gauss equation which can be written schematically: 

R^R + A'^. 

Thus, we obtain: 

[A, curl]B = RdB + dAdB + Ad'^B + A^dB. 

This yields: 

E = -(-A)-^(RaS + dAdB + Ad'^B + ^^95) + {-A)-^{AdA + (6.14) 
Using the elhptic estimate (6.7) on E^, we have: 
\\d'^E\\ 3, , < ||A£;|| 3, , 



3| 



+ U'^9B\\ 3, ^ + \\AdA\\ , . + \\A'\\ 



% ||R'l|Lf'L2(St)||95||z,^ooj;,6(St) + ||-4||ioc^6(S()||9^-B||L^°OL2(St) 

+ (PI|l-L6(S0 + ||9S||i^i6(so)||5^||L-L2(SO 

+Plli»L6(s,)(l|5^llLr^'^(s*) + ll^llLri«(s*)) 

where we used in the last inequahty the Sobolev embedding (6.5) on E^, the bootstrap 
estimates (5.5) for A, the bootstrap estimate (5.3) for R and the estimates (6.10) for B. 
Together with the Sobolev embedding (6.4) on Ej, we finally obtain: 

||9£;|Uoo^3(s,) + \\d'E\\^^^.^^^^ < 

Next, we estimate ||ii^||^2j^oo(-Sj). We first claim the following non sharp embedding on 
Ef. For any scalar function f on E^, we have: 

||(-A)-^'u||Locfs,) < Ikll 14, . + 1^11 13, (6.15) 

IIV J y^t) II lli-g-(s^) II llL-g-(St) ^ ' 

The proof of (6. 15) requires the use of Littlewood-Paley projections on E^ and is postponed 
to Appendix A. We now come back to the estimate of ||-E'||l2^oo(S()- Using (6.14) and 
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(6.15), we have: 

\\E\\L'lL^(T.t) 

where we used in the last inequahty the bootstrap assumptions (5.5) for A, the boot- 
strap assumption (5.3) for R, the bootstrap Strichartz estimate for B,see (5.16), and 
Proposition 6.4. This concludes the proof of the lemma. □ 



7. Estimates for DB 

The goal of this section is to derive estimates for OB, with B = A~^curl (A) using 
the wave equation (3.32) satisfied by each component of Ai. We provide the proof of two 
important propositions concerning estimates for Hcurl A and HB, with B — A'^curl (A). 
The proofs makes use of the special structure of various bilinear expressions and thus is 
based not only on the bootstrap assumptions for Aq,A, k and R but also some of our 
bilinear bootstrap assumptions. 

We will need the following straightforward commutation lemma. 

Lemma 7.1. Let (p be a so{3, 1) scalar function on Ai. We have, schematically, 

djincf)) - D{dj(f)) = 2(A^),- dxd^cl) + doAQd(t> + A.^d(t>. (7.1) 

We also have: 

[□,A](/. = -4r^VaV6(9o0)+4n-^V„nV(.(ao(ao</.))-2Vor''V„V(.(/. (7.2) 

= dAo + A^, 
F(2) = ddAo + AdA + A^, 

where Va and denote induced covariant derivatives on Ej applied to the scalars 4>, do4> 
and 9o(9o0). 

Remark 7.2. The derivation of (7.2) involves the full use of the Einstein equations. The 
resulting structure of the terms on the right hand side of (7.2) is crucial to the strategy of 
the proof of the main result and reflects "hidden" null cancellations. We refer the reader 
to (7.27) -(7.34), where this structure allows us to use bilinear estimates. 
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Proof. We start with the following general covariant calculation for any scalar function (p 
on M: 

[D^,q0 = O. (7.3) 
This follows trivially from the vanishing of the spacetime Ricci curvature, i.e. 

On the other hand, Lemma 3.4 yields: 

(e,)m(D0)^ = □(9,0) - 2(A^), B^d^cf) - D^(A,), ^9^0 - (A;,), ^(A^)^ '^9,0. 

Together with our Coulomb like gauge condition, wc obtain for djcf), j = 1, 2, 3: 

^,(□0) - □(a,-0) = 2(A^),- dxd^(l> + do{Aoyj 9^0 + A^dcj), 

which proves the first part of the lemma. The proof of the second part of Lemma 7.1 is 
postponed to Appendix C. □ 

7.1. Estimates for Dcurl{A). 

Proposition 7.3. The following estimate holds true, 



^||(-A)-in(mr/(A),)|U^(A4) < MV. 



i=l 



Proof. We have: 

n{djiA,) - d,{A,)) = d,{n{A,)) - diiD{A,)) + [D,d,]{A,) - [□,9,](A,). (7.4) 

We evaluate the first term on the right-hand side of (7.4) by differentiating (3.32). We 
obtain: 

dj{D{Ai)) = -d,{di{doAo)) - dM'diAi) + a,(/if )), (7.5) 



where h\^^ is given by: 



/if ^ = A^diAi + AodA + ASAq + 



We estimate h\^^ using the bootstrap assumptions (5.5) and (5.6) for A and ^40, the 
Sobolev embedding on (6.9), and the Bochner inequality (6.8) on S^: 

||(-A)-ia,(/.f))|U.(^) < \\h^^\\LHM) + \\Ah^''\\,^j^^^^ (7.6) 



< 



\\A^dl{Ai)\\L2(^_M) + ||^o||L2ioo(S,)||9A||iooi2(s< 

+||^IU?°l6(s,)||9^o||l°°l3(S4) + ||A|||ooi6(s,) 
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In view of (7.5), wc have: 

a,(n(A,))-5,(n(A,)) = -d,{d,{dM) + d,{d,{d^Ao)) (7.7) 

-d,{A'd^Ai) + di{A'djAi) + d^{hf^) - di{hf) 

— 'Hj ' 

where hf^ is given by: 

hf^ = Q,,{A\ ^i) + ^5(5oAo) + A^dK + a,(/if )) - d,{hf), 
and where the quadratic form Q^j is defined as Qij{(f),ip) = di(j)djip — djCpdic/). Note that 
the most dangerous term in h^f is QijlA^Ai). Using the bihnear assumption (5.14), the 
Sobolev embeddings on (6.9) and (6.5), the bootstrap assumptions (5.5) and (5.6) for 
A and Aq, and the estimate (7.6), we have: 

||(-A)-^(/^g))|U.(^) (7.8) 

;-A)-|(g,,(A\ A,))|U.(^) + \\AdidoA,)\\^,,j^^^^ + \\A'dA\\^,^.^^^^ 



+ ||(-A)-^9,(/^f))|U.(^) + \\{-A)-W,{hf)h.^M) 
< MV + ||A|U^.^6(s^)||a(aoAo)||^^^^3^^^^ + ||A||ioo^e(s,)||aA|Uc«^.(s,) + M^e^ 

Next, we consider the commutator terms in the right-hand side of (7.4). In view of 
(7.1), we have: 

[□, dj]{A,) = 2{A^)j ^ d,d,{A,) + (7.9) 

where h^^^ is given by: 

hfj^ =doAodiAi) + A^d{Ai). 

Using the Sobolev embeddings on (6.9) and (6.5), and the bootstrap assumptions (5.5) 
and (5.6) for A and ^40, we have: 



< \\(_A)-k(h^^^ 



^-Hf^lhHM) (7.10) 



< II^AII 

< Mh'. 

Next, we consider the term {A^)j'^ dxd^{Ai). We have: 

(A^)j^ dxd^(Ai) = -(Ao)j^ dodi(Ai) + (Ao)jododo(Ai) (7.11) 
+ {A')j ^ didrniAi) - {A')jo didoiA) + A^dA. 
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Note that the most dangerous terms in (7.11) are the third and the fourth one. They will 
both require the use of bilinear estimates. 

We deal with each term in the right-hand side of (7.11), starting with the first one. We 
have: 

= di{AodA) + di{Ao)dA + A^dA, 

which together with the Sobolev embeddings on (6.9) and (6.5), and the bootstrap 
assumptions (5.5) and (5.6) for A and Aq yields: 

\\{-A)-^{Ao),^ dommLHM) (7.12) 
< IIAaAIU.(^) + \mAo)dA\\^^j^^^^ + ||A^9A||^^^^e^^^^ 

^ ll^0||L2ioo(S4)||9A||iooi2(s,) + ||9i(Ao)||Looi3(St)||9A||iooi2(s,) 

+ ll^llL^L6(St)II^A|U?°^^(St) 

Next, we consider the second term in the right-hand side of (7.11). For that term, we 
would hke to factorize the do derivative in order to get two terms of the type do{Aodo{A)) 
and do{Ao)do{A), and then conclude using elliptic estimates and Sobolev embeddings on 
Ej. A similar strategy worked for the first term in the right-hand side of (7.11). But it 
does not work directly for this term since (— A)~29o is not necessarily bounded on L^(Et). 
Thus, we first start by showing how one may replace one do with d. Using the identity 
(3.3) relating A and R, we have: 

{Ao)jo dodo{A,) = iAo),odo{do{Ai)) + A^dA 

- {Ao),o do{d,{Ao)) + AodoiRor-.) + A.'^dA 
= Aod,{do{Ao)) + AoDoRo... + A2R + A^QA. 

Using the Bianchi identities for R, we have: 

DoRoi-- — D;R;j. .. 

Thus we obtain: 



{Ao)jo dodo{Ai) = Aodi{do{Ao))+AoBiRi,.. + A''R + A''dA 

= Aodi{do{Ao)) + diiAoR) + di{Ao)R + A'R + A'dA. 
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Using the Sobolev embeddings on Et (6.9) and (6.5), the bootstrap assumptions (5.5) and 
(5.6) for A and ^40, and the bootstrap assumption (5.3) for R, we obtain: 

||(-A)-^((Ao),o <9o<9o(A,))|U2(^) (7.13) 

< \\Aod,{do{Aom^,j^^^^ + \\AoR\\lhm) + \mAo)R\\^^^.^^^^ 
+ ||A2r|| , + IIA^^AII 6, , 

< \\M\L^Leij:,)\\di{do{Ao))\\^^^3^^^^ + \\Ao\\l2loo(j:^)\\R\\l^l^j:,) 

+ ||9i(^o)||L-L3(St)l|R'IU?°^^(St) + ll^llz,f>L6(Et)(II^IU?°^'(St) + II^A||iooi2(s,)) 

Next, we consider the third term in the right-hand side of (7.11). We have: 

(A')/™ a;a^(A) = A%XdiiA)) + A^dA 

= diA^diiA)) + d{A^)diiA) + A^dA. 

Together with the bihnear assumptions (5.9) and (5.15), the Sobolev embeddings on 
(6.9) and (6.5), and the bootstrap assumptions (5.5) and (5.6) for A and Aq, we obtain: 

||(-A)-^((A'),-9,9^(A,))|U.(^) (7.14) 

< \\A%{A)\\,.^M) + \\{-A)-Hd{A^)di{A))U.^M) + IIA^^AII^^^e^^^^ 

< Mh'. 

Finally, we consider the fourth term in the right-hand side of (7.11). We would like 
to factorize the Oq derivative in order to get two terms of the type do{A''di{A)) and 
dQ{A^)di{A), and then conclude using the bilinear assumptions (5.9) and (5.15). A similar 
strategy worked for the third term in the right-hand side of (7.11). But it does not work 
directly for this term since {—A)~2do is not necessarily bounded on L^(Et). Thus, as for 
the second term, we first start by showing how one may replace do with d. Using the 
identity (3.3) relating A and R, we have schematically: 

ao(A)-5i(^) + A2 = Ro,.. 

which yields: 

{A')jodidoiA,) = {A')jodi{do{A,)) + A'dA 

= iA')jo diidiiAo)) + iA')jo di{Roi..) + A'dA 

= di{{A')jo Roi..) + di{A')R + Ad''Ao + A''dA 

= di{{A^)jo Roi.-) + A.^R + A.''dA + Ad^Ao, 

where we used in the last inequality our Coulomb like gauge choice which yields di{A'') — 
V/(A') = A^. Thus, we have: 

(A')jo dido(Ai) = di((A')jo Roi-) + hf, (7.15) 
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where h\f is given by: 



hlf = A^R + A^^A + Ad^Ao. 



Using the Sobolev embeddings on (6.9) and (6.5), the bootstrap assumptions (5.5) and 
(5.6) for A and Aq, and the bootstrap assumption (5.3) for R, we obtain: 

\\{-A)-Hh^^)U.^M) (7.16) 
< IIA^RII 6, +\\A^dA\\ 6, ^ + \\Ad^Ao\\ 6, , 



Next, we estimate the first term on the right-hand side of (7.15). Since Roioo = 0, the 
terms Roj .. are of two types: Roimn or R-oiOm- Now, from the symmetries of R and the 
Einstein equations, we have: 

Also, in view of the hnk between R and A (3.3), we have schematically: 

R-mnOi — 9m{An) — dn{A^) + A and Rninm — ^n(^i) ~ di{A„) + A . 

Thus, we obtain schematically: 

Roi.. = dA + A^. 
which, using the Coulomb gauge, yields: 

di{{A')joRoi..) = A'did{A) +A^dA 

= d{A^di{A)) + d{A^)di{A) + A^dA. 

Together with the bilinear assumptions (5.9) and (5.15), the Sobolev embeddings on Ej 
(6.9) and (6.5), and the bootstrap assumptions (5.5) and (5.6) for A and Aq, we obtain: 

||(-A)-i(aK(^O.oRo^.))IU^(A.) (7.17) 

< \\A%{A)U.^M) + \\{-i^rHd{A^)di{A))\\^.^M) + ||A^aA||^^^6^^^^ 

< Mh\ 
Now, (7.15)-(7.17) imply: 

||(-A)-^((AO,o ^MA^))\WiM) < Mh\ (7.18) 
Finally, (7.11)-(7.18) imply: 

||(-A)-^((A^)/ d^d,{Ai))\W^M) < (7.19) 
In the end, (7.4), (7.7)-(7.10), and (7.19) yield: 

||(-A)-ln(9,(A))||L^(Ai)<MV. 
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This implies: 

\\{-^)-hn{d,{A,)-dUMWiM)<Mh\ 

which concludes the proof of the proposition. □ 

7.2. Estimates for \3B. Here we derive a wave equation for each component oi B = 
A~^curl (A) and prove the following, 

Proposition 7.4 (Estimates for DB). The components Bi = {—A)^^{curl{A)i) verify 
the following estimate, 

3 

J2\\daBi\\L^M) < (7.20) 



1=1 



We also have, 



J2\\do{do{B.mL^^M) < Me. (7.21) 



Proof. The estimates for OB are simpler than those for dOB and OdB. We prove first 
the estimates for dOB and derive those for OdB using the commutation formula (7.1). 
We have: 

D{Bi) = [D,{-Ar']{curl{A),) + {-A)-\D{curl{A),)) 

= -{-A)-'[D,A]{-A)-\curl{A),) + {-A)-\n{curl{A)i)) 

= -{-A)-'[a,Am) + {-A)-\n{curl{A)i)). (7.22) 

Thus, we obtain: 

\\dn{B,)\\L2^M) (7.23) 

< ||a(-A)-ip,A](A)IU^(A.) + ||5(-A)-^(n(cKH(A),))IU^(^) 

< ||(-A)-^[n,A](s,)|U.(A,) + ||(-A)-^(n(ci.w(A)0)|U^(^) 

< \\i-A)-^n,Am)U.^M) + Mh', 

where we used Proposition 7.3 in the last inequality. 

In view of (7.23), we need to estimate ||(— A)~2[n, A](Sj)||i2(jvi)- Recall the commu- 
tator formula (7.2): 

[□, A](/) = -4k'''VaV„{do<j)) + 4n-^V6nV6(9o(9o0)) - 2VoA;«''V„V;,0 (7.24) 

= dAo + A^, 
= ddAo + AdA + A\ 
Using the bootstrap assumptions (5.5) for A and (5.6) for ^40, we have: 

||i^^'^||LrL3(E.) < WdAoh^LS^^,) + ||A||io.^e(s,) < Me, (7.25) 



THE BOUNDED CURVATURE CONJECTURE 49 

and: 

< WddAoW^^^s^^^^ + ||A||ic«i6(s,)||aA||iooi2(s,) + ||A||iooi6(s,) 

< Me. 

Using (7.24)-(7.26) together with the estimates of Proposition 6.4, we obtain: 

||(-A)-^P,A](50l|L^(>t) (7.27) 

< ||(-A)-ir''V„V,(ao(A))]|U^(>,) + \\{-A)-^n-'V,nV,{do{do{Bi)m\L^(M) 
+||(-A)-^[Vor^V„V6(S0]|U^(>,) + ||F«|UrL3(E,)||a'(S0IU2(A^) 

< ||(-A)-^[r^V,V5(9o(50)]||L^(Ai) + ||(-A)-5[n-iVfenV5(9o(9o(i?/)))]||L2(.M) 

+ ||(-A)-^[Vor^V,V,(50]|U2(>i) + MV + Me||5o(9o(fiO))llL2(A^). 

Next, we estimate the various terms in the right-hand side of (7.27). The first and 
the third will require bilinear estimates, while the second will require the estimate Vn e 
L°°{M). We start with the first one. We have: 

k'^'VaVbidoiBi)) = V4r^Vfe(9o(A))] - v„r'V6(ao(50) 

where we used the constraint equations (2.2) for k in the last equality. Together with the 
Bochner inequality on (6.8) and the bilinear assumption (5.10) , we obtain: 

\\{-A)-'^k^'WaV,ido{Bim\LHM) < \\k^%{do{Bi))]\\LHM) < (7.28) 
Next, we estimate the second term in the right-hand side of (7.27). We have: 

n-^Vi,nV,{do{do{Bi))) = V'[n-^V ,ndo{do{Bi))] - {n'^ An - n-^\Vn\^)d^{do{Bi)). 

Together with the estimates (6.1) for the lapse n and the Sobolev embedding on (6.9), 
this yields: 

\\{-A)-"^[n-^VmVb{do{do{Bi)))]\\mM) (7.29) 

< \\n-^Vmdo{d^{Bi))\\L^^M) + \\{n-^An - n-^\Vn\^)do{d^{Bi))\\^,^,^^^^ 

< (\\Vn\\Lo^ + \\n~^An - n~'^\Vn\^\\Lf^L^(Y.t))\\dQ{do{Bi))\\L2(^j^-) 



< 



Me\\doidoiBi))\\L2^My 



Remctrk 7.5. Note that there is no room in the estimate (7.29). In particular, the esti- 
mate II Vn||ic=o(^) < Me given by (6.1) is crucial as emphasized in remark 6.1. 
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Finally, we consider the third term in the right-hand side of (7.27). Recall from (2.5) 
that the second fundamental form satisfies the following equation: 

Vokab = E,, + Fif, (7.30) 
where E is the 2-tensor on defined as: 

Eab = R-aObO, 

and where F^^^ is given by: 

In view of the estimates (5.7) for k and (6.1) for n, F^^^^ satisfies the estimate: 

lli^if l|LrL3(E0 < W'^Ml LrLHj:.) + ml^LH^,) < Me. (7.31) 

Next, we consider the term involving E in the right-hand side of (7.30). Using the maximal 
foliation assumption, the Bianchi identities and the symmetries of R, we obtain: 

V^-Eafe = D"Ra06o + AR = — D°Roo6o + AR = — (9o(Roo6o) + AR = AR 

which together with the bootstrap assumptions (5.5) for A and (5.6) for Aq, and the 
bootstrap assumption (5.3) for R yields: 

W^'^^'^'^krLh^.) ^ ^ (7.32) 

Now, we have: 

E,,VaVb{Bi) = V"K5V,(50] - V"E,bVb(50 

which together with the bilinear estimate (5.11), the estimates of Lemma 6.4 for B and 
(7.32) yields: 

\\{-A)-^EabVaVb{Bi)]\\L2^M) (7-33) 

< \\{-A)-'^V^[E,,V,(Bi)]\\l.^m) + \\(-A)-^V^E,,d,(Bi)]\\r^.^M) 

< \\'RaObodb{Bi)\\L2(^j^) + \\V"-Eab\\^^^3^^^^\\d{Bi)\\LooL6(^St) 

< Mh\ 

(7.30), (7.31), (7.33) and the estimates of Lemma 6.4 for B yield: 

||(-A)-|[Vor''V,V„(i?z)]|U2(^) (7.34) 

< \\{-A)-'^[E,,VaVb{Bi)]\Wm + ||(-A)-^[F„?)V„V,(A)]IU^(X) 

< M'e' + \\F^^\\LrLH^.)\\d'm\LrL^ii:.) 
< 

Finally, (7.27)-(7.29) and (7.34) yield: 

||(-A)-^p,A](S,)|| <MV + Me||9o(9o(S,)))|| 
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which together with (7.23) imphes: 

\\dU{B,)\\L2^M) < + Me\\do{do{Bi)))\\L^(M)- (7-35) 

Recalhng (3.25), we have: 

doidoiB,)) = -n{B.i} + A{B.i) + n-Wn ■ V(5,), 

which together with the estimates of Lemma 6.4 for B, the estimates (6.1) for n, and 
(7.35) yields: 

WdoidoiBmLHM) (7.36) 

< \\n{Bi)\\LHM) + ||A(5,)||l2(^) + \\Vn-V{Bi)\\LHM) 

< + Me\\do{dom))\\LHM) + \\d\B,)\\L^LHi:,) + \\Vn\\L^\\d{B,)\\L^LHi:,) 

< + Me\\doidoiBm\LHM)- 

Choosing e > such that Me is small enough to absorb the term ||(?o((9o(i?i)))||L2(_vj) in 
the right-hand side, (7.35) and (7.36) gives the desired estimate for both ||c?ni?||i:,2(_yvi) 
and ||9o(9o(-Bj)))||i2(^) of the lemma. 

□ 

8. Energy estimate for the wave equation on a curved background with 

bounded curvature 

Recall that Cq = T, the future unit normal to the foliation. Let tt be the deformation 
tensor of cq, that is the symmetric 2-tensor on A4 defined as: 

In view of the definition of the second fundamental form k and the lapse n, we have: 

TTab = -2A;a6, TTaO = TToa = n'^VaU, TToo = 0. (8.1) 

In what follows T-L denotes an arbitrary weakly regular null hypersurface^^ with future 
normal L verifying g{L,T) — —1. We denote by ^ the induced connection on the 2- 

surfaces "H fl S^. 

We have the following energy estimate for the scalar wave equation: 

Lemma 8.1. Let F a scalar function on M., and let 0o and <pi two scalar functions on 
Eq. Let 4> the solution of the following wave equation on A4: 

0|so = 00, <9o0|s„ =01- 

Then, (p satisfies the following energy estimate: 

||^0||L-L2(St) + SUp(||y0||i2(^) + ||L(0)||i2(^)) 

n 

< ||V</)o|U2(Eo) + ||0l|U2(Eo) + ||i^lUw, (8.3) 



32= 



i.e. it satisfies assumptions (5.1) and (5.2) 



52 



SERGIU KLAINERMAN, IGOR RODNIANSKI, AND JEREMIE SZEFTEL 



where the supremum is taken over all null hypersurfaces H satisfying assumptions (5.1) 
and (5.2). 

Proof. We introduce the energy momentum tensor Q^js on Ai given by: 

In view of the equation (8.2) satisfied by (f), we have: 
Now, we form the 1-tensor P: 

Pa — QaOi 

and we obtain: 

where vr is the deformation tensor of cq. Integrating a specifically chosen region of A4, 
bounded by Eg, and we obtain: 



< 



< 



l|V0o||i2(So) + ||0l||i2(So) + 



Fdo(t)dM 



M 



+ 



l|V0o||i2(Eo) + II0i|Il2(Eo) + II^IU2(X)||«9o0||l2(A4) 



M 



Next, we deal with the last term in the right-hand side of (8.4). In view of (8.1), we 
have: 



M 



= -2 / QabnM + / n-'V'nQoidM 
'm 



-2 f da<j)db(pk''''dM + I iigk {g^'^df^cpd^cj)) dM + [ n-^V''nda(i)do<j)dM 
Jm J m J 



= -2 / da(t>d^(t>k''''dM + 
Jm 

where we used in the last inequality the maximal foliation assumption. Together with the 
bilinear bootstrap assumption (5. 12) and the estimates (6.1) for the lapse n, this yields: 



JM 



< 



\\ka-d''(l)\\L^M)\\d(l)\\LHM) + l|Vn||ioo(^)||a(/)||i2(^) 



H 



which together with (8.4) concludes the proof of the lemma. 



□ 
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Remark 8.2. The most dangerous term in the right-hand side of the previous inequality 
is \\ka.d^(l)\\L'^(M)- Usually, when deriving energy estimates for the wave equation, this 
term is typically estimated by: 

||^a.9"0||L2(M) < \\k\\L^L^{i:,)\\d(f)\\Lt°L^{i:t) 

which requires a Strichartz estimate for k. This Strichartz estimate fails under the as- 
sumptions of Theorem 2.2, and we need to rely instead on the bilinear estimate (5.12). 

We have the following higher order energy estimate for the scalar wave equation: 

Lemma 8.3. Let F a scalar function on A4, and let 0o o-nd 0i two scalar functions on 
Eq. Let (p the solution of the wave equation (8.2) on M.. Then, satisfies the following 
energy estimate: 

11^(^0) lUrL^CE.) + ||5o(ao0)|U2(^) +SUp {mdct>)U^H) + ||L(a0)|U2(^)) 

■H 

< ||VVo||l2(So) + ||V(/.i|U2(So) + \\VF\\l2^m), (8.5) 

where the supremum is taken over all null hypersurfaces T-L satisfying assumption (5.1) 
and (5.2). Furthermore, □(9,0) satisfies the following estimate: 

IP(9,-0)|U2(^) < Me(||V^0o|U^(Eo) + l|V0i|U2(So) + ||9F|U.(^)). 

Proof. We derive an equation for Differentiating (8.2), we obtain: 

r n{d,<j,) = d,F + [n,d,M, 

\ 9i0|so = 9j0o, 9o(9j0)|so = {dj{do(l)) + [do, dj](l)) |so = dj(j)i + A0i + AV^q. ^ ' ' 
Applying the energy estimate of Lemma 8.1 to (8.6), we obtain: 

\\9{dj(f))\\L^L^j:,) + sup{\\f{dj(f))\\L2(n) + \\L{dj(p)\\L2(H)) 

< II V(a,0o)IU2(So) + \\djcf>i + A01 + AV0o||l2(So) + \\djF + [□, OMl^m)- 
which after taking the supremum over j = 1, 2, 3 yields: 

11^(90) IUrL^(EO + sup (11^(90) |U.(«) + ||L(a0)|U2(^)) (8.7) 

H 

^ ||VVo|U2(Eo) + I|V0i||l2(So) + ||VF||i2(_yK) + ||A0l||i2(So) + 1 1 A V0O 1 1 L2 (So) 

+ \\^9(t>\\L^L^Et) + SUpllp, 9j]0||l2(^), 

i 

where the term Ad(j) in the last inequality comes from the commutator formula (3.22) 
applied to [do, dj]. 

Next, we estimate the last term in the right-hand side of (8.7). In view of the commu- 
tator formula (7.1), we have: 

[n,dj](f) = 2{A%^ dxd,,(f) + doAod(f) + A^d(f) (8.8) 
= A'di{di<f>) + A%{docf>) + A°9o(9,0) + (^°),o9o(9o0) + doAod<f> + A^a^ 
= A'di{di<l>)+A%{do(f>) + hj, 
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where hj is defined in view of tfie identity (3.17) as: 

hj = A'^do{di(l)) + n-^Vjudoidocp) + doAodcj) + A^dcj). 

We estimate the various terms in the right-hand of (8.8) starting with hj. The Sobolev 
embedding on (6.5), the bootstrap estimates (5.6) for Aq and (5.5) for A, and the 
estimate (6.1) for the lapse n yield: 

Whjh^M) (8.9) 

^ ll^°l|L?L-(E,)||9o(9i0)||L-L2(EO 

+||Vn||Loo(^)||9o(9o0)||L2(_yK) + ||9o^o|U-l3(s,)||90||l-l6(s,) + ||A||ioci6(S4)||90|U-L6(St) 
< MeiWdoidomLHM) + ||9(a0)|UrL^(E.))- 

Note again in view of the previous inequality that the estimate Vn G L°°{A4) is crucial 
as emphasized by Remarks 6.1 and 7.5. Next, we deal with the first and the second term 
in the right-hand of (8.8). Using the bihnear estimate (5.13), we have: 



\\A%{dicj>)\\mM) + \\A%{docl>)\\mM) < Me l^sMWfidMLHH) + lim</')IU^(H) j , 
which together with (8.8) and (8.9) yields: 

\m,dj]<f>\\L^M) (8.10) 

< Me (sup{\\f{di<P)\\L2^n) + mdomLHu)) + Me{\\do{domLHM) + \\d{d<l>)h^r.^j,,)). 



H 

It remains to estimate the term ||y(9o0)||L2CH)- Let us define the vectorfield N = L — cq. 
Since g(L, cq) = —1, and since L is null, is tangent to E^. Decomposing eo — L — N, 
we obtain schematically: 

iWo0)i < mvN<jy)\ + mm)\ (8.11) 

|y(iV,40)| + |9(L(0))| 

\f{dcl>)\ + |L(a0)| + |(D7V)(a</-)| + |(DL)(a0)| + \Adcl>\ 

which together with the assumptions (5.1) and (5.2) for H, and the embedding (10.2) on 
H yields: 

11^(^00) (8.12) 

< mamLHu) + \\L{dmLHH) + mN){dmLHH) + mmamLHu) + \\Ad<ph.^n) 



< 

< 



< mdcl>)\\mn) + \\L{dcf>)U.^n) + (l|D7V|U3(«) + ||DL|U3(^) + \\AUs^n))\mLHn) 



< 
< 



(II W)|U.(H) + ||L(90)|U.(«))(1 + ||VA|Uc«^.(s,)) 
||y(90)|U.(^) + ||L(90)|U.(«), 
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where we used the bootstrap assumptions (5.5) for A in the last inequahty. Finally, 
(8.10)-(8.12) yield: 



sup||[n,a,](0)|U2(^) < Me(^sup(||y(a0)|u + ||L(a0)|U2(^)j (8.13) 

+Me{\\do{domL-^iM) + ||5(90)|U,^L2(so)- 

Now, (8.7) and (8.13) imply: 

119(90) lUoo^.(s^) + sup (11^(90) |U.(«) + ||L(a0)|U2(^)) 

H 

^ ||V^0o||l2(So) + ||V0i||l2(2„) + \\VF\\l2(m) + ||A0i||l2(2„) + ||AV0o||l2(Eo) 

+ ||Aa0|Uoc^2(s,) + Me (^sM\\f{dcl>)\\LHn) + 11^(90) |U.(«)) 

+Me(||9o(9o0)|U2(^) + ||9(a0)|Uco^.(s,)), 

which together with the Sobolev embedding on E^, the bootstrap estimates (5.6) for Aq 
and (5.5) for A, and the fact that we may choose e such that Me is small enough, yields: 

||5(90)IUrL2(so +sup(||y(a0)|U2(H) + \\Lid<P)\\L^n)) (8.14) 

< I|v20o||l2(so) + l|V0i|U2(s„) + + Me||ao(ao0)|U2(^). 

In view of (8.14), we need an estimate for 9o(9o0). Proceeding as in (7.36), we have: 

L^(M) ^ IP(0)IU2(A4) + ||A0||i2(_vi) + ||Vn ■ V(0)||l2(_a4) (8.15) 

^ \\F\\l^M) + l|V^0||L-L2(St) + ||Vn||Lcx=||V0||L-L2(s,) 

^ II^IIl2(a^) + ||9V||l,-l2(S()- 

Finally, (8.6) and (8. 13)- (8. 15) yields: 

\\n{djcl>)\\L^M) < \\djF\\L^M) + mdML.^M) 

< \\dF\\L2^M) + Me (snpimdmn + 11^(90) |U2(^)) 
+Me(||9o(9o0)|U2(^) + ||9(a0)|Uo.^.(s^)) 

< Me(||V^0o|U^(Eo) + l|V0i|U2(s„) + ||9F|U2(^)) 

which together with (8.14) and (8.15) concludes the proof of the lemma. □ 

9. Proof of Proposition 5.7 

Here we derive estimates for R, Aq and A and thus improve the basic bootstrap as- 
sumptions (5.3), (5.4), (5.5) and (5.6). 
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9.1. Curvature estimates. We derive the curvature estimates using the Bel- Robinson 
tensor, 



Let 

Then, we have: 



Qafj-yS — 

Ra 7'^R/3A5(T+ *Ra 7'^*R/3A5(t 



Pa — Qa/37(5Co^O^O- 



D'^P^ = 3Q 

where tt is the deformation tensor of cq. We introduce the Riemannian metric, 

hap = gap + 2(eo)a(eo)/3 
and use it to define the following space-time norm for tensors U : 



(9.1) 
(9.2) 



\U\ 



Uai—akUa[—a'i^h ^ i 



Given two space-time tensors U, V we denote by U ■ V a given contraction between the 
two tensors and by \U -Vl the norm of the contraction according to the above definition. 

Let be a weakly regular null hypersurface with future normal L such that g{L, T) — 
— 1. Integrating (9.1) over a space-time region, bounded by So,St and Ji, and using 
well-known properties of the Bel-Robinson tensor, we have: 



|RP + 



H 



|R • ^ l|R|li2(Eo) + 



JM Jm 



We need to estimate the term in the right-hand side of the previous inequality. Note that 
since ttoo = 0, n^j = n'^VjU, and TTij — kij, the bootstrap assumption (5.3) for R, and 
the estimates (6.1) for n yield: 



/ w+ / 

Jn 



|R-^I' < e^ + llVnl 



|R'lli-L2(St) + 



M 



< e' + {Mef 



M 



The last term on the right-hand side of the previous inequality is dangerous. Schemat- 
ically it has the form | kYi? | . Typically this term is estimated by: 



kB? 



M 



< 



Z,2Loo(S()||R||Looi2(s,), 



requiring a Strichartz estimate for k which is false even in flat space. It is for this reason 
that we need the trilinear bootstrap assumption (5.19). Using it we derive. 



IR-LI 



(9.3) 



which, for small e, improves the bootstrap assumptions (5.3) and (5.4). 
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9.2. Improvement of the bootstrap assumption for Aq. Recall (3.31): 

AAo = AdA + A-\ (9.4) 

Using the elliptic estimate (6.7) and the Sobolev embedding (6.4) together with (9.4), we 
have: 

\\dAohrLS^^,) + \\d'Ao\\^^^s^^^^ < IIA^II^^^^s^^^^ (9.5) 

where we used the bootstrap assumptions (5.5) on A and the Sobolev embedding (6.5) in 
the last inequality. 

Next, using the Sobolev embedding (6.15) together with (9.4), we have: 

I|A||l|loc(e,) = ||(-A)-i(A9A + A=^)|U.^oc(s,) (9.6) 
^ l|A9A||^,^.^^^^ + ||Ai^,^.^^^^ 

< Mh\ 

where we used the bootstrap assumptions (5.5) on A and the Sobolev embedding (6.5) in 
the last inequality. 

Next, we consider ^o^o- In view of (9.4), we have: 

A(9oAo) = 9o(A9A) + ao(A=^) + [ao,A](^) 

= 9oAaA + a(AaA) + A^aA + a](Ao). 

Together with (C.4), we obtain: 

A(ao/lo) = 9oA9A + 9(AaA) + A2aA+[ao,A](/lo) (9.7) 
= /1 + 5/2, 

where /i is given by: 

/i = ^oA^A + A^SA, 

and where /2 is given by: 

h = AdA. 

In view of the bootstrap assumptions (5.5) for A and (5.6) for Aq, we have: 

||/i|U?°Li(s*) + ll/2|lic.ii(s,) (9-8) 

We will use the following elliptic estimate on S^: 
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Lemma 9.1. Let v a scalar function on St satisfying the following Laplace equation: 

Av^f, + df2. 
Then, we have the following estimate: 

lklU3(s,) + ~ ll/lll^HS*) + ll/2|li§(s,)- 

The proof of Lemma 9.1 requires the use of Littlewood-Paley projections on and is 
postponed to Appendix B. We now come back to the estimate of OqAq. In view of (9.7), 
Lemma 9.1 and the estimate (9.8), we have: 

WdoAohrLHi:,) + \\9doAo\\^^^^^^^^ < M\\ (9.9) 

Finally, (9.5), (9.6) and (9.9) lead to an improvement of the bootstrap assumption (5.6) 
for ^0- 

9.3. Improvement of the bootstrap assumption for A. Using the estimates for \3Bi 
derived in Lemma 7.4, the estimates for B on the initial slice Eq obtained in Lemma 4.3, 
and the energy estimate (8.5) derived in Lemma 8.3, we have: 

||5'5|U^i2(2,) <e + MV. (9.10) 
Using (9.10) with Lemma 6.5, we obtain: 

WdAWL^LHn,) < Wd'BU^LH^^) + WdEU^L.^^^) < e + (9.11) 
Next, we estimate Oq^A). Recall that: 

9o(A,) = 9,(A)+Ro,-... 

Thus, we have: 

||c?0^||L°°Z,2(St) % ||c?^o||L°°L2(Et) + ||R-||z,°°L2(St), 

which together with the improved estimates for R and Aq yields: 

||5o^||lc«l2(s,) <6 + (Me)i. (9.12) 

Finally, (9.11) and (9.12) lead to an improvement of the bootstrap assumption (5.5) for 
A. 

Finally, (9.3), (9.5), (9.6), (9.9), (9.11) and (9.12) yield the improved estimates (5.20), 
(5.21), (5.23) and (5.22). This concludes the proof of Proposition 5.7. 
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10. PARAMETRIX for THE WAVE EQUATION 

Let u± be two families of scalar functions defined on the space-time A4 and indexed 
by a; e S^, satisfying the eikonal equation g"^daU±dfsu± — for each a; e We also 
denote ^u±{t, x) — u±{t, x, u). Wc have the freedom of choosing '^u± on the initial slice 
So, and in order for the results in [42] , [44] to apply, we need to initialize '^u± on Sq as 
in [41]. The dependence of '^u± on u is manifested in particular through the requirement 
that on So the behavior of u± asymptotically approaches that oi x ■ uj. 

Let % i^u^ denote the null level hypersurfaces of '^u±. Let '*'L± be their null normals, 
fixed by the condition g( '*'-L±, T) = ^1. Let the vectorfield tangent to Sj "^N^ be defined 
such as to satisfy: 

= ±eo + "iV±. 

We pick ( ^e±)A-i A = 1,2 vectorfields in S^ such that together with ^N± we obtain 
an orthonormal basis of S^. Finally, we denote by y± derivatives in the directions 
(-e±)A, ^ = 1,2. 

Remcirk 10.1. Note that from the results in [43] (see Theorem 2.15 and section 3.4 in 
that paper) 'H<^u± satisfy assumptions (5.1) and (5.2). 

We record the following Sobolev embedding/trace type inequality on "H^ for functions 
defined on M., derived in [43] (see sections 3.5 in that paper). 

Lemma 10.2 (An embedding on % [43]). For any null hypersurface Hu, defined as above, 
and for any T^-tangent tensor F, we have: 

mWm < IIVF||^^^^3(^^^ + \\F\\LrLH^.), (10.1) 

and for any 2 < p < 4: 

\\F\\lp{h) ^ ||VF||^oo^2(s^) + ||F||ioc=^2(j]^). (10-2) 
For any pair of functions f± on IR^, we define the following scalar function on M.: 

p poo p poo 

^[f+J.]{t,x)= / e'^""+(*'^V+(Ac^)A'rfAcia;+ / / e'^ /-(Ac^)A'dAda;. 
We appeal to the following result from [42] [44] : 

Theorem 10.3 (Theorem 2.11 in [42] and Theorem 2.17 in [44]). Let 0o and 0i two 
scalar functions on Sq. Then, there is a unique pair of functions {f+,f-) such that: 

V'[/+,/-]|so = 00 anc?9o(V' [/+,/-]) I So = 0i- 
Furthermore, f± satisfy the following estimates: 

||A/+||l,2(R3) + ||A/_||i2(R3) < ||V0o||l2(So) + ||0l|U2(So), 

and: 

||A^/+||l2(R3) + ||A^/_||i2(]K3) < ||V^0o||l2(So) + ||V0l||i2(So). 
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Finally, □^[/+,/_] satisfies the following estimates: 

||n^[/+,/_]|U2(^) < Me(||V0o||L^(Eo) + WMlH^o)): 

and: 

||an^[/+,/_]|U.(^) < Me(||V^0o||L^(So) + l|V0i|U2(So)). 

Remark 10.4. The content of Theorem 10.3 is a deep statement about existence of a 
generalized Fourier transform and its inverse on Eq, and existence and accuracy of a 
parametrix for the scalar wave equation on M., with merely L"^ curvature hounds assump- 
tions on the ambient geometry. The existence of f± and the first two estimates of Theorem 
10.3 are proved in [42], while the last two estimates in Theorem 10.3 are proved in [44]. 

We associate to any pair of functions 0i on Eq the function ^om[0O) 0i] defined for 
(t, x) e M as: 

^omi^O^l] =V' [/+,/-] 
where (/+, /_) is defined in view of Theorem 10.3 as the unique pair of functions associated 
to (00) 0i)- In particular, we obtain: 

I|A/+||l2(m3) + ||A/_||l2(k3) < ||V0o||l2(So) + ||0i||l2(So), 

||A^/+|U2{M3) + ||A^/_||l2(]r3) < ||V^0o||l2{So) + ||V0i||l2(So), 

||OTom[0o,0i]||L2(^) < Me(||V0o|U2(So) + ||0i|U2(So)), (10.3) 

and: 

||9OTo„[(/)o,0i]||l-(ai) < Me(||V'0o||L2(Eo) + W^ML^i^,))- (10.4) 

Next, let ^'*M± two families, indexed by a; e and s e M^^, of scalar functions on 
the space-time M. satisfying the eikonal equation for each a; e §^ and s e M. We have 
the freedom of choosing on the slice E^, and in order for the results in [42] [44] to 

apply, we need to initialize '^''^u± on E^ as in [41] . Note that the families ^u± correspond 
to '^'^M with the choice s = 0. For any pair of functions /-t on M^, and for any s e M, we 
define the following scalar function on M.: 

^s[f+J-]{t,x,s)^ / e'^"''''+^''''^f+{Xuj)X^dXduj+ / e'^"'°"-(*'^V-(Aa;)A2ciAcia;. 

J§2 Jo J§2 Jo 

We have the following straightforward corollary of Theorem 10.3: 

CoroUciry 10.5. Let s e R. Let 0o ond (pi two scalar functions on E^. Then, there is a 
unique pair of functions (/+,/_) such that: 

'^s[f+, /-]|e, = 00 and 9o(V's[/+, /-])|e, = 0i- 
Furthermore, f± satisfy the following estimates: 

I|A/+||l2(r3) + ||A/-||l2(jj3) < ||V(/>o||l2(s^) + ||(/'i||l2(s,), 

33ln fact, we need to restrict s to an open interval of [0, 1] for which the bootstrap assumptions (5.3) 
(5.4) on R hold 



THE BOUNDED CURVATURE CONJECTURE 61 

and: 

Finally, □■^^ [/+,/_] satisfies the following estimates: 

\\D^s[UJ-]\\mM) < Me(||V0o||L^(E.) + ||0l|U^(E.)), 

and: 

\\dnM+,n\\L^(M) < Me(||V^0o|U^(E.) + ||V0i|U2(s,)). 

Next, for any s e R, we associate to any function F on the function ^{t, s)F defined 
for {t, x) e M as: 

M/(t,s)F = ^,[/+,/_](t) 

where (/+, /_) is defined in view of CoroUary 10.5 as the unique pair of functions associated 
to the choice {(po, (pi) = (0, —uF). In particular, we obtain in view of Corollary 10.5 and 
the control of the lapse n given by (6.1): 

||A/+||l2(k3) + ||A/_||l2(k3) < ||F||i2(s,), 

||A^/+||l2(R3) + ||A^/-||l2(ir3) < ||VF||i2(s,), 

s)FU.^M) < M6||F|U2(s,), (10.5) 

and: 

\\dn^{t, s)F\\l2^m) < Me|| VF||i2(s.). (10.6) 
Remctrk 10.6. Note that we have 

□ (^^ ^{t,s)F{s)ds^ ^F{t) + n^{t,s)F{s)ds. 

Now, we are in position to construct a parametrix for the wave equation (8.2). 

Theorem 10.7 (Representation formula). Let F a scalar function on M., and let 0o o-nd 
4>i two scalar functions on Sq. Let 4> the solution of the wave equation (8.2) on M.. Then, 
there is a sequence of scalar functions {(f)^^\ F^^^) , j >0 on M., defined according to: 

0(°) = *on^[0O,0l]+ t^{t,s)F^'\s,.)ds, FW=F 

Jo 

and for all j > 1: 

(f)U) = j ^(t, s)F^^\s, .)ds, F(^) = -D^^^-^) + F(^-i) 
Jo 

such that, 

and such that (f)^^ and F^^^ satisfy the following estimates: 

W^^Hl^l^i:,) + WF^^^Wl^m) < (Me)^(||V0o||L^(Eo) + II0iIU^(Eo) + \\F\\lhm)), 
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and: 

Proof. Let us define: 

F(°) = F and = *„^[0o, 0i] + /* 

Jo 

Then, we define iteratively for j > 1: 

F(^) = -□^(^■-i) + and 0^^) = [ ^(i, 

In view of Remark 10.6, note that for j > I: 

= □ / ^{t, s)F^^\s, .)ds = F(^) + /" .)ds, 
Jo Jo 

which yields: 

FU+i) = _f n^(t,s)F^^\s,.)ds. 
Jo 

Thus, we obtain in view of (10.5) and (10.6): 

<Me\\F^^^U^M), 

and: 

\\dF^'^'^\\,^.^M)<Me\\dF^^^\\mM)- 
Therefore, we obtain for all j > 2: 

\\F^'%HM)<iMey-'\\F^'%HM), (10-7) 



and: 



Also, we have: 



This yields: 



WdF^^^hHM) < {Mey-'\\dF^'^\\LHM)- (10.8) 

□0^°^ = F^°^ + n*om[</'o,0i]+ [ n^{t,s)F{s,.)ds, 

Jo 



F« = -n^om[(l>0, 0l] - r □*(^, S)F{S, .)ds 

Jo 

which together with (10.3), (10.4), (10.5) and (10.6) implies: 

\\F^'^\lhm) < Me(||V0o|U2(Eo) + WMlh^o) + \\F\\v^iM)), 

and: 

\\dF^'^\\L^(M) < MedlWolU^CEo) + l|V0i|U2(s„) + PFW^^^m))- 
Together with (10.7) and (10.8), we obtain for any j > 1: 

WF^'^LHM) < (Me)^'(||V0o|U2(Eo) + l|0ilU2(Eo) + \\F\\lhm)), (10-9) 
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and: 

WdF^'^hHM) < (M6y(||VVo|U^(Eo) + W^^MlHjio) + \m\LHM))- (10.10) 

We now estimate ^^■'^ j > 1. For j > 1, ^^■'^ satisfies the following wave equation: 

□0(i) ^ pU) _ pU+i)^ 



0(^)|eo = O, 9o(0(^)|eo = O. 
which together with Lemma 8.3, (10.9) and (10.10) yields: 

W^'^hrL^i^*) ^ iMey{\\vM\LH^o) + WMlh^o) + WnL^M)), (lo.n) 

and: 

\\d{dcl>(^^)\\LrL^i:,) < (M6)^d|VVo||L^(Eo) + l|V0i|U2(s„) + WdFU^M)). (10.12) 
Now, we have: 

□ ( = X](F(^) - F(^+i)) = F - 

\i=o / j=o 

which together with (10.10) and (10.12) yields in the limit j — >■ +oo: 
Note also that 

0(°)|so = 0oandao0(°)|Eo=0i, 

while for all j > 1, we have: 

0(^)|so = Oandao0(^)|so = O. 
Thus, Yl'j^ 0^''^ satisfies the wave equation (8.2), and by uniqueness, we have: 

+ 00 

This concludes the proof of the theorem. □ 



11. Proof of Proposition 5.8 (part 1) 

The goal of this and next sections is to prove Proposition 5.8. This requires the use 
of the representation formula of Theorem 10.7. In this section we derive the improved 
bilinear estimate (5.24), (5.25), (5.26), (5.27) and (5.28) of Proposition 5.8. We also derive 
the improved trilinear estimate (5.33). 
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11.1. Improvement of the bilinear bootstrap assumptions I. We prove the bihnear 
estimates (5.24), (5.25), (5.26), (5.27), (5.28). These bihnear estimates aU involve the 
L'^{M.) norm of quantities of the type: 

C(C/,90), 

where C(C/, dcf)) denotes a contraction with respect to one index between a tensor U and 
d(f), for a solution of the scalar wave equation (8.2) with F, 0o and 0i satisfying the 
estimate: 

||VVo||l2(Eo) + l|V0i|U2(So) + \\dF\\L2^M) < Me. 
In particular, we may use the parametrix constructed in Lemma 10.7 for 0: 

+ 0O 

j=0 

with: ^ 

0^°^ = *om[0O, 0l] + / *(i, S)F{S, .)ds, 

Jo 

and for all j > 1: 

(f)(^) = I '^{t,s)F^^\s,.)ds. 
Jo 

Thus, we need to estimate the norm in L^(A^) of contractions of quantities of the type: 

+ 00 „t 

C(t/,a(vl/„^[0o,0i])) + J] / C{U,d{^{t,s)F^^\s,.)))ds. 

j=o 

After using the definition of and ^{t,s), and the estimates for F^^^ provided by 
Lemma 10.7, this reduces to estimating: 

/ / C(t/,5(e^^""+(*'")))/+(Aa;)A'(iA(ia;+ / / C:(f/, 9(e^^ '^"-(*'^)))/_(Aa;)A'dAda;, 

where f± in view of Theorem 10.3 and the estimates for F, 4>q and 4>i satisfies: 

||AV±I|l2(m3) <Me. 

Since both half-wave parametrices are estimated in the same way, the bilinear estimates 
(5.9), (5.10), (5.11), (5.12) and (5.13) all estimate the norm in L^(A^) of contractions of 
quantities of the type: 



/poo 
,2 Jo 



where / satisfies: 

l|AVlU2(M3)<Me. (11.1) 

Now, we have: 
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and the gradient of on is given by: 

where '^b — |V( ^u)\~^ is the null lapse, and 

V ""u 
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IV ""ul 



is the unit normal to H •^u^^t along E^. Thus, the bilinear estimates (5.9), (5.10), (5.11), 
(5.12) and (5.13) all reduce to L^(Al)-estimates of expressions of the form: 



C[U, /] := / / e'^ ""b-^CiU, ''N)f{Xuj)X^dXduj, 
JS2 Jo 

where / satisfies (11.1). 

To estimate C[U, f] we follow the strategy of [21]. 



(11.2) 



< 



< 



ix -«(t,a.)/(Aa;)A^rfA 

r r+oo 

7§2 Jo 



(11.3) 

duo 



[ sup| 








I sup 









^Uj)\\L2duJ 



where we used Plancherel in A and Cauchy Schwarz in u. Now, since '^u has been 
initialized on Eq as in [41], and satisfies the eikonal equation on M, the results in [43] 
(see section 4.8 in that paper) under the assumption of Theorem 2.10 imply: 

sup II '^6"^||l-(>1) < 1- 
Together with the fact that / satisfies (11.1), and with (11.3), we finally obtain: 



.2 Jo 



< 



(11.4) 



L2(X) 



It remains to estimate the right-hand side of (11-4) for the contractions appearing in 
the bilinear estimates (5.24), (5.25), (5.26), (5.27) and (5.28). Since all the estimates in 
the proof are uniform in a;, we drop the index uj to ease the notations. 
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Remark 11.1. In the proof of bilinear estimates (5.24), (5.25), (5.26), (5.27) and (5.28), 
the tensor U appearing in the expression C{U, N) is either R or derivatives of solutions 
of a a scalar wave equation. In view of the bootstrap assumption (5.4) for the curvature 
flux, as well as the energy estimate for the wave equation in Lemma 8.1, we can control 
||C(t/, A^)||x,cx.L2(-^^) as long as we can show that C{U,N) can be expressed in terms of, 

Yl- L,f(j) and L{(j)). 

In other words, our goal is to check that the term C{U, N) does not involve the dangerous 
terms of the type: 

a and Lcj) 

where L is the vectorfield defined as L = 2T — L, and a is the two tensor on n Hu 
defined as: 

11.1.1. Proof of (5.24). Since A = curl (B) + E in view of Lemma 6.5, we have: 

\\A^dM)\\LHM) < ||(ctxr/(5))^9,(A)|U2(^) + ||E|U.^.(s^)||M|U^^^.(s,) (11.5) 
< \\icurliB)yd,{A)\\LHM)+Mh^ 

where we used in the last incquahty Lemma 6.5 for E, and the bootstrap assumption 
(5.5) for A. Next, we estimate \\{curl {B)y dj{A)\\L2(^j\^y Recall that we have: 

{curl {B)ydj{A) ^ejmn dm{Bn)dj{A). 

We are now ready to apply the representation theorem 10.7 to B. Indeed, according to 
Lemma 7.4, and proposition 6.4, we have 

DB = F, WdFU^^M) < (11.6) 

l|as(o)|u.(so) + ||a2s(o)|U2(s„) + ||a(aoS(o))|U2(s„) < Me. 

We are thus in a position to apply the reduction discussed in the subsection above and 
reduce our desired bilinear estimate to an estimate for, 

C{U,N) = ejm.N^dj{A) 

Now, we decompose dj on the orthonormal frame A^, /a, A = 1, 2 of E^, where we recall 
that fA,A — 1,2 denotes an orthonormal basis of Hu n E^. We have schematically: 

+ % (11.7) 

where ^ denotes derivatives which are tangent to Hu H T,f Thus, we have: 

C{U, N) =G,„. N^N,dN{A) + f{A) = f{A), 

where we have used the antisymmetry of Ejm- in the last equality. Therefore, we obtain 
in this case: 

\\mN)\\L^r.Hn.)<mA)\\L^L^^n.y 
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It remains to estimate ||y(A)||iooi2(-^^). Since we liave A = curl {B) + E in view of Lemma 
6.5, we obtain: 

\\f{A)\\L^LHH^) < \\f{dB)\\L^LHn^) + \\fiE)\\L^L^^n^) 

< mdB)h^L.^n.) + l|5i?IUj»L3(E0 + \\d'E\\^^^s^^^^ 

< mdB)U^^.^n.) + Me, 

< Me 

wliere we used tlie embedding (10.1) and tlie estimates for E given by Lemma 6.5. Fur- 
thermore, we have in view of Proposition 7.4 and Lemma 8.3 the following estimate for 
B: 

liy(as)|u.^.(^„)<Me. 

We finally obtain: 

||y(^)IU5fL2(«„) <Me. 
This improves the bilinear estimate (5.9). 

11.1.2. Proof of (5.25). In view of Lemma 6.5 A — curl (B) + E. Arguing as in (11.5), 
we reduce the proof to the estimate of: 

\\{curlByd,{dB)\\L2^M)- 
Since B satisfies the wave equation (11.6), the quantity C{U,N) is in this case, 

C{U,N) =e,^.N^d,{dB). 
Using the decomposition of dj (11.7) and the antisymmetry of Ejm-, we have schematically: 

ejm-Nmdj{dB) = Ejn,. N^NjdN{dB) + f{dB) (11.8) 

= f{dB) 

= f{dB)+fidoB) 

= f{dB)+L{dB) + {'DL + 'DN + A)d{B), 

where in the last equality we used the decomposition (8.11) for y(9o(-B)). Together with 
the assumptions (5.1) and (5.2) on Hu, and the Sobolev embedding (10.2) on Hu, we 
obtain: 

< mdB)\\L^n.) + \\LidB)\\L^n.) 

+ (I|D^IU3(H„) + I|D^I|l3(H„) + P|U3(^„))||a5|U6(^„) 

< (1 + WAWl^lhj:.) + \\dAh^L^^Mf{dB)U2^n.) + ||L(5i?)|U2(^„)) 

< mdB)h^H.) + \\L{dB)\\L^n.), 

where we used the bootstrap assumption (5.5) for A in the last inequality. Now, we have 
in view of Lemma 7.4 and Lemma 8.3 the following estimate for B: 

\\f(dB)U^r.^n^) + \\L{dB)U^r.HH.) < Me, 
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which improves the bihnear estimate (5.10). 

11.1.3. Proof of (5.26). Since B satisfies a wave equation in view of Lemma 7.4, the 
quantity C{U,N) is in this case: 

iVjRoj.. — RoAT..- 

Thus, using the fact that L = T + N, L = T — N and the symmetries of R, we deduce: 

NjRoj.. = -Rll-- 

which together with the bootstrap assumption for the curvature flux (5.4) improves the 
bihnear estimate (5.11). 

11.1.4. Proof of (5.27). We have kj. = and A = curl (B) + Em view of Lemma 6.5. 
Arguing as in (11.5), we reduce the proof to the estimate of: 

\\{curlBydj(j)\\L2(M)- 

Since B satisfles a wave equation in view of Lemma 7.4, the quantity C{U, N) is in this 
case: 

Using the decomposition (11.7) for dj and the antisymmetry of &jm-, we obtain schemat- 
ically: 

which improves the bilinear estimate (5.12). 

11.1.5. Proof of (5.28). We have A = curl (B) + E in view of Lemma 6.5. Arguing as in 
(11.5), we reduce the proof to the estimate of: 

II (cur/ Bydj(j)\\L2^M)- 

Since B satisfies a wave equation in view of Lemma 7.4, the quantity C{U,N) is in this 
case: 

Using again the decomposition (11.7) for dj and the antisymmetry of &jm-, we obtain 
schematically: 

which improves the bilinear estimate (5.13). 
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11.2. Improvement of the trilinear estimate. In this section, we shall derive the 
improved trilinear estimate (5.33). Let Qa/SjS the Bell- Robinson tensor of R: 

We need an trilinear estimate for the following quantity 



M 



We have A — curl (B) + £■ by Lemma 6.5. Arguing as in (11.5), we reduce the proof to 
the estimate of: 

Q.j^s{curl{B))je2el 



M 



Making use of the wave equation (11.6) for S we argue as in the beginning of section 11.1 
to reduce the proof to an estimate of the following: 



^iX '^u{t,x) u)y- 



'M Jq 

where / satisfies: 



-jm- 



""NmQj...) f{\uj)\^d\dujdM 



||AVlU2(M3)<Me. 



Arguing as in (11.3) (11-4), we obtain: 

JM Js'^ Jo 

[ ^ (ejm. ""NrnQj ■ ■ ■) ( f^^ ""(*'^)/(Aa;) A^dA 

7S2 \Jq 



< 



< 



< 



J§2 



duo 



du 



sup II ""6 Il-(AI) ) ( SUp IIGj^. ""NrnQj-.-WL^^^L^Hu.^) 



\X^f{\u)\\Lidu: 



where we used Plancherel in A and Cauchy Schwarz in ui. Thus, we finally obtain: 



M 



3,3 



< supllej™. Ar^g,-...||i2^ L\'H.^)Me + M-'e 



(11.10) 



Next, we estimate the right-hand side of (11.10). Since all the estimates in the proof 
will be uniform in cj, we drop the index u to ease the notations. The formula for the 
Bell-Robinson tensor Q yields: 

Qj... = R j • -R-A-- + dual 

— — — RjL..R. L.. — L-R-L-- + RjA-R-A-- + dual^ 
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where we used the frame L, L, fA,^ = 1,2 in the last equahty. Thus, we have schemati- 
cally: 

&jm- NmQj--- — R(R • L+ Ejm- NmRjA-) 

Decomposing ej with respect to the orthonormal frame N, fB,B — 1,2, we note that: 

On the other hand, decomposing Rba-- further and using the symmetries of R, one easily 
checks that Rs^.. must contain at least one L so that it is of the type R • L. Thus, we 
have schematically: 

ejm-NmQj...^^{^-L). (11.11) 

Thus, in view of (11.10), making use of the bootstrap assumptions (5.3) on R and (5.4) 
on the curvature flux, we deduce. 



M 



< (Me)' + Me\\RRL\\LiLHH.) 

< {Mef + Me||R|U2(^)||Ri||i5.i2(^„) 



JM 



In other words, 

^-^-^ < {Mef. (11.12) 
which yields the desired improvement of the trilinear estimate (5.19). 

12. Proof of Proposition 5.8, (part 2) 

In this section we prove the bilinear estimates II. We start with a discussion of the 
sharp Strichartz estimate. 

12.1. The sharp L^(A^) Strichartz estimate. To a function / on and a family "^u 
indexed by a; e of scalar functions on the space-time M. satisfying the eikonal equation 
for each a; e we associate a half- wave parametrix: 

Js^ Jo 

Let p be an integer p and ip he a smooth cut-off function on supported on the interval 
[1/2,2]. We call a half- wave wave parametrix localized at frequencies of size A ~ 2^ the 
following Fourier integral operator: 

/ / e'^ ""(*'^V(2"^A)/(Aa;)A2(iAda;. 

i§2 Jo 

We have the following L^(A^) Strichartz estimates localized in frequency for a half wave 
parametrix which are proved in [45] : 
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Proposition 12.1 (Corollary 2.8 in [45]). Let f be a function on M.^, let p E N, and 
let ip be as defined above. Let be a family of scalar functions on the space-time M. 

satisfying the eikonal equation for each lo G and initialized on the initial slice Eq as in 
[41] . Define a scalar function 4>p on Ai as the following oscillatory integral: 

(t>p{t,x)^ / / e^^ ""(*'^V(2~^A)/(Aa;)A'dAda;. 

Then, we have the following L^[M.) Strichartz estimates for (pp.- 

UpWlhm) < 2§||V'(2-^A)/||^.(M3), (12.1) 

WdMLHM) < 2^||V'(2-^A)/||^2(j,3), (12.2) 

Wd'MLHM) < 2t||^(2-fA)/||i2(i,3). (12.3) 

Note that these Strichartz estimates are sharp. 

12.2. Improvement of the non sharp Strichartz estimates. In this section, we 
derive the improved non sharp Strichartz estimates (5.31) and (5.32). In order to do this, 
we first estimate the LjL''(St) norm of dB using the L*{Ai) Strichartz estimate together 
with Sobolev embeddings on E^. 

Corollciry 12.2. B satisfies the following Strichartz estimate: 

WdBh.^r^,,,) < Me. 
Proof. Decompose B as before, with the help of Theorem 10.7, 



+ 00 

Thus is suffices to prove for all j > 0: 

W^L^mj:,)<{Mey^\ (12.5) 

The estimates in (12.5) arc analogous for all j. so it suffices to prove (12.5) in the case 
j = 0. In view of the definition of the estimates for B on the initial slice So obtained 
in Lemma 4.3, the estimate (7.20) for dOB, and the definition of "^om and ^'(t, s), (12.5) 
reduces to the following estimate for a half wave parametrix: 

< ||AV||l^(m3). (12.6) 

Next, we introduce (p and ip two smooth compactly supported functions on M+ such 
that ip is supported away from and: 

(p{X) + = 1 all A e M. (12.7) 

p>0 



d 



f{\u)\^d\dLO 
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We define the family of scalar functions (pp for p > — 1 on as: 

/ / e^^ ""(*'^V(A)/(Aa;)A'dAdc^, (12.8) 

and for all p > 0: 

(t>p{t,x)= / e^^""^*'^V(2"^A)/(Aa;)A'cZAda;. (12.9) 
In view of (12.7), we have: 



which yields: 



The estimate for 0_i is easier, so we focus on (f)p for p > 0. Using the Sobolev embedding 
(6.6) on Ej, the L^(Al) Strichartz locahzed in frequency of Proposition 12.1, and the fact 
that is supported in (0, +oo), we have: 

4 3 
||50pI|l2L7(S«) < ll^</'plll4(M)ll'90plll4ioo(s,) 

- 9 - 

< (2^||^(2-^A)/|U.(M3))' (2f ||^(2-^A)/|U. 

< 2-^||A2V^(2-PA)/|U.(M3) 

< 2-&||AV||l2(m3). 
Together with (12.10), we obtain: 

d(f r e'^ "<''^^f{Xu)X'dXdw] <(j2 I II AVIU2(M3) < || AVIU^Cms), 

which is (12.6). This concludes the proof of the Corollary. □ 

Lemma 6.5 and Corollary 12.2 yield: 

ll>l|U?L^(EO < I|55||l?L^(E,) + ll^l|L?Lr(E,) < M6 + (12.11) 

which is an improvement on the bootstrap assumption (5.16). 
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12.3. Improvement of the bilinear bootstrap assumptions II. In this section, we 
derive the improved bihnear estimate (5.29) and (5.30) of Proposition 5.8. Recall the 

decomposition A = curl (B) + i? of Lemma 6.5. Using the bootstrap assumption 5.5 for 
A, the estimates for E given by Lemma 6.5 and the Sobolev embedding on (6.5), we 
have: 

\\i-ArHdAdE)h^M) + \\{-Ar^dEdE)U^M) 

Together with the decomposition A = curl (B) + E oi Lemma 6.5, this implies that the 
proof of the bilinear estimates (5.14) and (5.15) reduces to: 

||(-A)-i(g,,(ciiW(S),ciiW(S)))|U.(^)+||(-A)-^(9(ciiW(S)')9KciiH(S)))|U2(^) <M^ 
where the bilinear form Qij is given by Qij{(l),ip) = di(()djip — djCpdiip. Also, note that: 

d{curl {By)di{curl (B)) = Q,j{dB,dB) + AdB d^B. 

Together with the Sobolev embedding (4.7) on Ej, the bootstrap assumptions (5.5) on A, 
and the estimates of Lemma 6.4 for B, we obtain: 

-A)-^d{curliBy)diicurl m)\\LHM) 

:-A)-kQ,,{dB, dB))h2^M) + II (- A)-^ {A dB d'B) U^^m) 

-A)-^Q,^{dB,dB))U.^M) + \\AdBd^B\\^^^.^^^^ 

;-A)-^(Q,,-(9S,9S))||i2(^) + ||A||z,,^i6(s,)||ai?|Uj»L6(E,)||a2S|U^^L2(s,) 
;-A)-5(Q^^.(a5,a5))|U2(^) + MV. 
Finally, the proof of the bilinear estimates (5.14) and (5.15) reduces to: 

\\{-A)-HQ,,{dB,dB))\y^M) < (12.12) 

Next, we focus on proving (12.12). Decomposing B according to Theorem 10.7, we 
have: 

+ 00 

||(-A)-^(g,,(as,as))|U.(>,) < \\{-^)-HQ^A^^^'^\^4>^^^))\\L^M)■ (12.13) 

m,n=0 

Thus it suffices to prove for all m, n > 0: 

||(-A)-i(g,,(a</.(-),90(^)))|U.(^,) < {Mer+\MeT+\ (12.14) 

The estimates in (12.14) are analogous for all m, n, so it suffices to prove (12.14) in the 
case (m, n) = (0,0). In view of the definition of the estimates for B on the initial 



< 
< 



< 

< 
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slice So obtained in Lemma 4.3, estimate (7.20) for d\I\B, (12.14) reduces to the following 
bilinear estimate for half- wave parametrices: 



A)-ig,,(a0«,a(/.(^)) ^ < ||AVi|U^(R3)||AV: 



2 L2 



with, 



0W = 



S2 ^0 



fk{XLj)X^dXdu, k = l,2. 



(12.15) 



(12.16) 



Recall the smooth cut off functions (p and introduced in the proof of Corollary 12.2. 
We define two families of scalar functions 0^, j = 1, 2 for p > — 1 on as: 



p poo 

and for all p > 0: 

Js^ Jo 

In view of (12.7), we have: 



(12.17) 



(12.18) 



which yields: 



-A)-^g,,(a0«, 



,(1) ^^(2)^ 



L2 



< Yl \\i-^)-HQ^M\9<Pf^))h^. (12.19) 

The estimates involving (p^'^l are easier, so we focus on (p^^ , 0^^^ for p,q > 0. We may 
assume q > p. Note that the structure of Qij implies: 



which yields: 



< 
< 

< 
< 



||(-A)-^(g,,(a0«,a0f)))|U.(^) 
\\{-A)"^d{d'4^-d<p^^^)\\mM) + ||(-A)-i(^-aVi') 
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where we used the bootstrap assumption (5.5) for A in the last inequahty. Together with 
the L^{M.) frequency locahzed Strichartz estimate of Proposition 12.1, and the fact that 
■0 is supported in (0, +oo), we obtain: 

||(-A)-i(g,,(a0«,a0f))|u.(^) < 2^+t||^(2-^A)/i|U.(M3)||^(2-«A)/2|U.(M3) 

< 2i-5||A2^(2-PA)/i|U2(K3)||A2^(2-«A)/2|U2(M3). 

Since we assume q > p, this yields: 

J2 ||(-A)-|(Q,,(90«,90f ))|U.(^) (12.20) 



< 



p,g>-i 



J2 2-^||A2V(2-^A)/i|U2(m3)||AV(2-^A)/2|U2(m3) 



< $^||AV(2-^A)/i||i.(^3) $^||A^^(2-^A)/2||i. 
\p>-i / \?>-i 

^ ||A^/i||l2(i[53)||A^/2||l2(IR3). 

Finally, (12.19) and (12.20) imply (12.15). This concludes the proof of the improved 
bilinear estimates (5.29) and (5.30). 

Finally, (12.11), the results in section 11.1 and section 12.3, and (11.12) yield the 
improved estimates (5.24), (5.25), (5.26), (5.27), (5.28), (5.29), (5.30), (5.31), and (5.33). 
This concludes the proof of Proposition 5.8. 

13. Propagation of regularity 

The goal of this section is to prove Proposition 5.9. Recall from the statement of Propo- 
sition 5.9 that we assume that the estimates corresponding to all bootstrap assumptions 
of the section 5.3 hold for < i < T* with a universal constant M. For the convenience 
of the reader, we recall some of these estimates below 

||R'||Lj°°L2(St) + IIR- • L\\L2{n) ^ ^: (13.1) 

||>1||l-l2(eo + ||a>l|U-L2(s*) < e, (13.2) 

and: 

MolUf L2(St) + ll^^olUf L2(St) + ||^o||L2L°°(St) + II || L^L3(St) 

+ \\ddAo\\ 3, < e. (13.3) 

Note also that we have (6.1) with a universal constant M, i.e. 

\\n - + ||Vn||Loo(^) < e. (13.4) 
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Finally, for any weakly regular null hypersurface "H and any smooth scalar function (j) 
on M, 

\\kj-d^(l)\\L2(M) ^ esup||y0||i2(<^), (13.5) 
n 

and 

\\A^dj<f>\\L^M) < esup||y0|U2(^), (13.6) 
n 

where the supremum is taken over all null hyp ersurf aces T-L. 

We now turn to the proof of Proposition 5.9. The estimates will be similar to the one 
derived in Proposition 5.7 and 5.8. Thus, we will simply sketch the arguments without 
going into details. Note that the crucial point is to check that the null structure - which 
is at the core of the proof of Proposition 5.7 and 5.8 - is preserved after differentiation. 
This will be done in section 13.3. 

13.1. Elliptic estimates. Recall that we need to control ||DR||2,^i2(s^). Note that we 
may restrict our attention to the control of D^Rj^Q,^ and DoRijo,^. Indeed, all the other 
components are recovered from the Bianchi identities and the standard symmetries of R. 
Furthermore, since the estimates are for D^Rj^Q,^ and DoRjjQ/j similar, we will restrict 
only to the control of D^Rj^Q,^. These components, according to the Cartan formalism, 
can be expressed in the form 

Rija/J = {diAj - djAi + [Ai, Aj])^^ . (13.7) 

Since 

D^R = diR + AR 

we can estimate 

||D;Rija/3||L2(S() < ||(?iRija/3||L2(St) + ||R||L6(St) || A||i3(s,) . 

Furthermore, from (13.7) 

Proceeding the same way with all other components of DR, we derive 

I|5R||l2(so< l|5^A|U2(so+e'. (13.8) 

In view of (13.8), it suffices to derive bounds for ddA. Furthermore, since Aq 
satisfies an elliptic equation, and hence better estimates, we focus on the estimates for 
ddA. We sketch below the estimates just for d'^A since the estimates ddoA are similar. 
To estabfish bounds on ||5^A||i2(Sj) we use the following analog of Lemma 6.5: 

Lemma 13.1. The following decomposition 

OA = curl (dB) + E' 
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holds with E' satisfying: 

Furthermore, A satisfies 

Proof. Recall from Lemma 6.5 that we have the following decomposition 

A = curl B + E, 

with 

E = -{-A)-\RdB + dAdB + Ad'^B + A^dB) + {-/^)-^{AdA + A^) (13.9) 
given in (6.14). We now introduce the new variables dA and dB, hnked by the equation 

dA = curl (dB) + E', 

where 

E'^dE + [d, curl]B, 

so that 

E' = -d{-A)-^ {RdB + dAdB + Ad'^B + A^dB) + AdA + A^) + AdB. 
It then follows that 

and 

||5-E'||l?°L3(S,) < (^||-R|U-L6(S,) + ||9A||iooi6(S,) + ||>l||i-Ll2(S,)) ||5-B||iooi6(S,) 

+ ||^|U?°L6(Et)ll^^-6|U°°L6(St) 

+ ||^|Uf>Z,6(Et)||9^||L-L6(St) + ||^|li-L9(E0 

< e {\\dR\\L^L^i:,) + ||9M|U^^i2(s,) + ||9-^5|Uc»^2(2,)) + 

< e (||9i?'|U^^L2(E0 + ||9^i?||LrL^(so) + 

Using that in our, localized, setting the L^{^t) norm dominates the LF'^^t) norm we obtain 
that 

||9M|Uoo^.(S0<||a^i?||L,-L^(E0+e' 

and 

||9-E'||l°°l3(s,) < e||9^-B||Looi2(s,) + e^. 
The other estimates are proved in the same. This concludes the proof of the lemma. □ 

In view of Lemma 13.1, it remains to estimate d^B. This will be done following the 
same circle of ideas used for S^B. 
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13.2. The wave equation for dB. Recall from (7.22) that we have, schematically: 

□5- {-A)-^[n, A]B + {-A)-^n{curl A). 
Using the commutation formula (7.1) we obtain 

ndB = F (13.10) 

where F is given by 

F = A^dedB + A^d^B + A^dB + doA^OB + 9(-A)-^p, A]B + d{-A)-^n{curl A). 

(13.11) 

Using the energy estimates for the wave equation derived in Lemma 8.3 and the proof 
of the non sharp Strichartz estimate in section 12.2, we obtain the following proposition. 

Proposition 13.2. We have the following estimates 

• Energy estimate 

£ :=sup(||y92S|U2(^)+||L92S|U2(^)) + ||a92S|Uooi2(s,) < ||a92S|U2(s„)+||9F|U2(^)+e. 

(13.12) 

• Non sharp Strichartz estimate 

S := < + \\dF\\L2^M) + e. (13.13) 

Furthermore, we have the following estimate for ddodoB 

\\ddodoB\\L2^M) < ||aa2S|U2(So) + \\dF\\L^M) + e. (13.14) 

In view of Proposition 13.2, we need to estimate dF. This is done in the following 
proposition: 

Proposition 13.3. F satisfies the following estimate 

WdFWmM) < e {\\d'Bo\\mi:o) + ||5'i^i|U^(Eo) + S + S+ \\dF\\mM)) + e'. 

The proof of Proposition 13.3 is postponed to the next section. In view of Proposition 
13.2 and Proposition 13.3, we obtain the control for ||9i9^-B||L^L2(St) by its initial data 
which together with the elliptic estimates of the previous section yields the desired control 
for ||DR||ic,o^2(2t)- This concludes the proof of Proposition 5.9. 

13.3. Proof of Proposition 13.3. Prom (13.10), we have 

dF = dA^dedB + A^ddedB + d{A°d^B + A^dB + doA^dB) 
+d''{-Ay'[D,A]B + di'{-A)-'D{curlA). 
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This yields 

\\dF\\L^{M) < \\dA'd,dB\\L^M) + \\A'dd,dB\\L2(M) + \\P,^]B\\mM) (13.15) 

+ \\U{curl A)\\l2^m) + WdiA'^d^B + A^dB + doA^dB)\\L2^M) 

< \\dA%dB\\L2(^M) + \\A'dd,dB\\L2^M) + \\[D,A]B\\l2^m) 
+ \\n{curl A)\\l2^M) + l-o-t, 

where we neglect the cubic terms and the terms involving A^ since, as in the proof of 
Proposition 5.7 and 5.8, they are significantly easier to treat. 

Next, we isolate the terms D{curl A) and [□,A]i? on the right-hand side of (13.15). 
Prom the proof of Proposition 7.3, we have 

D{curl A) = d{A^diA) + Qij{A^, Ai) + terms involving Aq + cubic terms. 

We deduce 

\\n{curlA)\\L2^M) < \\A%dA\\L2^M) + \\dA%A\\L2^M) + \\Qij{^\ A,)\\l2^m) + l-o-t. 

(13.16) 

Next, we deal with [□, A]i?. According to (7.24), we have schematically 

[□, A]E = r^V,V5(9o5) + n-^V,n^,{do{d^B)) + VoA;'^'V„V65 
+ terms involving + cubic terms. 

We deduce 

||[n,A]S|U2(^) < \\k'^VaVb{doB)\\L2^M) + \\^^k'^yaVbB\\L2^M) 

+ \\n-^VbnVb{do{d^yB))\\L2^M) + l-O-t. 
Together with (13.15) and (13.16), we obtain 

\\dF\\L^M) < \\dA'dedB\\L2^M) + \\A'ddedB\\L2^M) + \\dA'deA\\L2^M) (13-17) 
+ \\A'd,dA)\\L2^M) + \UA^',Ae)\\mM) + l|A;"'V„V,(9o5)|U2(^) 
+ ||VoA;"''V„V6S||i2(^) + \\n-^VbnWb{do{doB))\\L2^M) + l-o-t. 
We will use the following bilinear estimates. 
Lemma 13.4. We have 

\\A'dd,dB\\L2i^M) + \\A'd,dA\\L2i^M) + ||A;"'V„Vfc(cio5)IU2(>i) 
< e(||a9^i?|U.(So)+^ + 5) + el 

Lemma 13.5. We have 

\\dA%dB\\L2^M) + \\dA'diA)\\L2^M) + \\QiM',A,)\\L2^M) + W^ok^'V oVbEWL^^M) 
< e{\\dd^B\y^^,) + ||9F|U.(^)) + e\ 
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The proof of Lemma 13.4 is postponed to section 13.3.1 and the proof of Lemma 13.5 
is postponed to section 13.3.2. We now conclude the proof of Proposition 13.3. In view 
of (13.17), Lemma 13.4 and Lemma 13.5, we have 

\\dF\\mM) < ||n-V6nV6(9o(9oS))|U2(^) + e(||a92S|U2(s„) + ||9F|U2(^)) + 
Using the estimates (13.4) and (13.14), we have 

\\n'^VbnVb{do{doB))\\L2(^M) ^ \\n'^Vn\\L--^M)\\ddodoB\\L2(M) 

< e(||a9^5|U.(s,) + ||aF|U.(^)) + e^ 

and we deduce 

\\dF\\L2^M) < e(||aa^S|U^(E„) +S + S+ ||9F|U2(^)) + 
which is the desired estimate. This concludes the proof of Proposition 13.3. 

13.3.1. Proof of Lemma 13.4- In view of the identity dA = curl {dB) + E' of Lemma 13.1, 
we have 

\\A'ddfdB\\L2^M) + \\A'dedA\\L2^M) + ||A:"'V„V,(ao5)|U2(^) 

< \\A'ddedB\\L2^M) + ||r^V„V,(9ofi)|U2(^) + /.o.t. 
Together with the bihncar estimates (13.5) and (13.6), we deduce 

\\A'ddedB\\L2^M) + \\A'dedA\\L2^M) + \\k'''VaVb{doB)\\L2^M) 

< €Sup\\fddB\\L2^n) + l-o.t. 

H 

Arguing as in (11.8), we finally obtain 

\\A'dd,dB\\L2^M) + \\A%dA\\L2^M) + ||^"'V„Vfe(ao5)|U2(^) 

< esnp{\\fd^B\\L2in) + \\Ld^B\\L2^n)) + l-o-t. 

n 

< e{\\dd'BU2(^^,)+£ + S) + e' 

which is the desired estimate. This concludes the proof of Lemma 13.4. 

13.3.2. Proof of Lemma 13.5. In view of the wave equation (13.10) satisfied by dB, we 
may use for dB the parametrix constructed in Lemma 10.7: 

+00 

dB = Y,cP^^\ 

j=0 



with: 



and for all j >1: 



0^°^=*om[</'O,0l]+ / ^{t,s)F{s,.)ds, 

Jo 

Jo 
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Furthermore (p^^^ and F^^^ satisfy the fohowing estimate: 

\\dd<l>^'^h^LHj:.) + WdF^^^h-^M) < 6^d|a^i^o||L^(Eo) + I|5'5i||l^(Eo) + WdFh^^M))- 

We will show that the proof of the bilinear estimates of Lemma 13.5 all involve the 
L^(A4) norm of quantities of the type: 

C{U,d{dB)), 

where C{U,d{dB)) denotes a contraction with respect to one index between a tensor U 
and d{dB). Using the parametrix for dB discussed above, and arguing as in section 11.1, 
we obtain the analog of (11.4): 

\\C{U,d{dB))\\L.^M) < {\\d'Bo\\mi:o) + \\d'B,\\,^.^^,) + \\dFU.^M)) (13.18) 



X (^sup||C(C/,7V)||^,(^) 

where the supremum is taken over all weakly regular null hypersurfaces, and where N is 
the unit normal to Pi inside E^. 

We are now ready to prove the bilinear estimates of Lemma 13.5. Using the decompo- 
sition of Lemma 13.1, we have 

\\dA%dB\\L2^M) + \\dA%A\\L2^M) + WQiM'^MlLHM) (13.19) 

< \\Qij{dB,dB)\\L2^M)+l-o.t. 

which is of the type C{U,d{dB)) with U — ddB. Now, arguing as in section 11.1, we 
have in this case 

C{U, N) ^Eij didBNj = fdB 
and we deduce from (13.18) and (13.19) 

\\dA%dB\\L2(^M) + \\dA%A)\\L2(^M) + \Uj{A\Ae)\\L^M) (13.20) 

< (^sup\\fdB\\L2^H)^ {S + S+ WFU^iM)) + 

< e{£ + S+\\F\\L2^M))+e''. 

Next, we consider the term Vak^-^V a^hB. Recall from (7.30) that we have 

'^okab = R-a060 + l-O.t. 

It follows that, 

L'^{M) ^ ||R'a06oVaV6-B||2,2(_;K) + l.0.t. (13.21) 

The right-hand side of (13.21) is of the type C{U, d{dB)) with U — Ra06o. Now, arguing 
as in section 11.1, we have in this case 

C(C/,A^) =Rjvo.. = R-iv 
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and we deduce from (13.18) and (13.19) 

\\dA'd,dB\\L2(M) + \\dA'dtA)\\L2(M) + \\QiM'^^i)\W{M) (13.22) 



|^sup||R • L||i2(^)^ {S + S+ ||-F||l2(^)) + 

< e{£ + S+\\F\\L2^M))+e' 



where we used (13.1) in the last inequahty. Finally, (13.20) and (13.22) yield the desired 
estimate. This concludes the proof of Lemma 13.5. 

Appendix A. Proof of (6.15) 

The goal of this appendix is to prove (6.15). We first introduce Littlewood-Paley 
projections on which will be used both for the proof of (6.15) and Lemma 9.1. These 
were constructed in [43] (see section 3.6 in that paper) using the heat flow on Ej. We 
recall below their main properties: 

Proposition A.l (Main properties of the LP Qj [43]). Let F a tensor on E^. The 
LP-projections Qj on E^ verify the following properties: 
i) Partition of unity 

J2Qj = i- (A.l) 

j 

a) W -houndedness For any 1 < p < oo, and any interval / C 

\\QiF\\l.(ji,) < (A.2) 
Hi) Finite band property For any 1 < p < oo. 



||Ag,F|U.(s,) < 22^-||F|U.(s,) 
||g,F||L.(E.) < 2-2^||AF|U.(s,). 



(A.3) 



In addition, the estimates 



hold together with the dual estimate 

||Q,-VF|U.(s,) < 2iF|U2(s,) 
iv) Bernstein inequality For any 2 < p < +oo and j & 1^ 

IIQ,i^lU.(E.) < 2^^'"'^'II^IIl^(e.) 

together with the dual estimates 
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We now rely on Proposition A.l to prove (6.15). Using Proposition A.l, we have for 
any scalar function v on Tit'. 



< J]2¥||g,(-A)-VIU^(E0 



< 



< 



J]2-i||g,/llL^(E, 



u>o / \j<o / 



13 , , 

z,-g-(St) 



- ii-^iIl^(e,) + ii-/i|l¥(so- 



This concludes the proof of (6.15). 

Appendix B. Proof of Lemma 9.1 
Recall that v is the solution of 

Av^h + df2. 

Using the harmonic coordinate system on of Lemma 6.2, we obtain: 

= {gij-5ij)d^^v + rdjv + rf2 + fi + d{f2) (B.l) 

= df^[{g,^ - s,j)v] + dj[rv] + dj{r)v + r/2 + A + ^(a), 

where d and A denote the derivatives and the fiat Laplacian relative to the coordinate 
system defined above, as opposed to the frame derivatives d and the Laplacian A, and 
where F is the corresponding Christoffel symbol. Now, we use the following standard 
elliptic estimates on M^: 

(-A)-^a, G C{L-2{W'),L\R^)), {-A)-% e £(Li(M'^)) 

(-A)-^ e jC{l\w'),l\r^)), (-A)-^4 e /:{l\r^),lI{r^)), 

where the notation C{X, Y) stands for the set of bounded linear operators from the space 
X to the space Y. Together with (B.l) and our assumptions on the harmonic coordinates 
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(6.2) (6.3), this yields in a coordinate Patch U: 

\MLHu) + \\dv\\j^^^ 

< Wl^U) + + \\vi\\LHU){\mLHU) + l|r||Le(C/)) 

+lk2||L4([;)(||9r||^4(^) + ||r||^^(j,)) + ||/i||li(c/) + \\M\^^^^^ 

where v = Vi + V2. We then sum the contributions of the covering of by harmonic 
coordinate patches U satisfying (6.2) (6.3). Eventually to increasing C{6), we obtain 



< <^(||^IU3(E,) + + C{S){\\v,\\lH^,) + ||^;2||l4(S*)) + ||/i||l1(S*) + ll/2||^i(s,)- 

Recall from Lemma 6.2 that we have the freedom of choice for 6 > 0. By choosing S > 
small enough, we obtain: 

lkl|L3(E0 + ^ C'(5)(||^i|U.(s,) + |k2|U4(E.)) + ||/i||l1(E.) + ll/2||^3(^^)- (B.2) 

5 > is now fixed. Thus, C{5) — C is a constant which may not be small. We now choose 

vi = Q>ov and V2 = Q<ov. 

Together with (B.2), we obtain 

||^^||l3(so + WMIl^^e,) ~ ll<9>o^lU^(St) + \\Q<ov\\l^(i:,) + ||/i|Ui(St) + ll/2||^i(s,)- ^^"^^ 

In view of (B.3), we still need to estimate ||(5>o'*^||L2(St) and ||(5<o'i^||L4(Et)- We start 
with ||(5>o'i^||L2(s^). We have: 

\\Q>ov\\mi:,) < ||(-A)-^Q>o/i|U2(E.) + ||(-A)-^Q>o9/2|U.(EO. (B.4) 

Next we estimate each term in the right-hand side of (B.4) starting with the first one. 
Using Proposition A.l, we have for any scalar function / on E^: 

||(-A)-ig>o/|Uoo(E.) < J2\\Q^(-^y'fh-(^^) 

j>0 

< 5^2¥||g^.(-A)-VlU^(Eo 

j>0 



< 



< 



< 



j>0 



E2- 
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Taking the dual, we obtain for fi. 

||(-A)-^Q>o/i|U2(s,) < ||/i|Ui(E,). (B.5) 

Next, we consider the second term in the right-hand side of (B.4). Using property i) of 
Proposition A.l, we have: 

||(-A)-^g>oaMU.(s,)< J2 m-^r'dQif2\\L^i^,y (B.6) 

j>o,iez 

We now estimate the right-hand side of (B.6). We consider the two cases j > I and j < I 
separately. If j > /, we obtain using Proposition A.l: 

\\Qj{-A)-'dQif2\\mi:,) < 2-'^\\dQif2\\mi:,) (B.7) 
If J < /, we obtain using Proposition A.l: 

||g,(-A)-iagj2|U^(Eo < 2-'||(-A)-i9(-A)^gj2lU^(so (b.s) 

< 2-'||(-A)-i9(-A)i|U(i.(so)||gj2||L^(E.) 

< 2-^||V^(-A)-i^(,.(^,))||/,||^3^^^^ 

where we used the fact that 

||V^(-A)-i£(i.(s,)) < 1 
thanks to the Bochner inequality on (6.8). Finally, (B.6)-(B.8) yields: 

||(-A)-^g<oa/,|U.(,,) < ( Y: 2-^^'^- + E 2-^) < ll/2||,i(^^). (B.9) 

\j>l,j>0 0<j<l / 

(B.4), (B.S) and (B.9) yield: 

||g>o^^lU^(E.) < ||/l|Ui(E,) + 11/211^3^^^^. (B.IO) 
In view of (B.3), we still need to estimate ||g<o'y||L4(Et)- We have: 

||g<oH|L4(E.) < ||(-A)-^g<o/i|U4(s,) + ||(-A)-^g<o9/2|U4(EO- (B.ll) 
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Next we estimate each term in the right-hand side of (B.ll) starting with the first one. 
Using Proposition A.l, we have for any scalar function / on E^: 

3J 



< 5^2^||g,(-A)-VIU^(E. 



j<0 



< J]2¥||Q,(-A)-V||, 



< 



< 



< 



'L§(St) 

j<0 

j<0 



L3(Et) 



L3(St)' 



Taking the dual, we obtain for fi. 



||(-A)-iQ<o/i||^4^^^^<||A|U.(s,). (B.12) 

Next, we consider the second term in the right-hand side of (B.ll). Using property i) 
of Proposition A.l, we have: 

||(-A)-^g<oa/2|U4(s,) < Yl \m-^r'dQif2\\LHi:,). (B.13) 

j<0,l<=Z 

We now estimate the right-hand side of (B.13). We consider the two cases j > I and j < I 
separately. If j > Z, we obtain using Proposition A.l: 

||g,(-A)-^agj2||L4(Eo < 2-t^||9gj2||L^(Eo (b.m) 

< 2-t^-+'||g,/2||L^(E.) 

< 2-2^+^ II fall 3 . 

If J < Z, we obtain using Proposition A.l: 

||g,-(-A)-i9gj2|U4(E0 < 2^2-'||(-A)-i9(-A)^gj2||L^(E.) (B.15) 

< 2^2-'||(-A)-^9(-A)i|U(^.(^,))||gj2|U^(E.) 



< 2¥-^||V^(-A)-i^(^.(^,))||/2||^3^^^^ 



3j I 



< 2^-2 II foil 3 , 

where we used the fact that 

||V2(-A)-i£(i.(s,)) < 1 
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thanks to the Bochner inequahty on Ht (6.8). FinaUy, (B.13)-(B.15) yields: 

||(-A)-g<oa/,|U4(,,) < [ J] 2-i^-+f + Y: 2"-^) 11/211,3^,^^ < ||/,||^3^^^^(B.16) 

\i<j<o j<i,j<o / 

(B.ll), (B.12) and (B.16) yield: 

||Q<0^IU4(E.) < ||/l|Ui(E.) + ll/2||^i(^^)- (B.17) 

Finally, (B.3), (B.IO) and (B.17) imply: 

This concludes the proof of Lemma 9.1. 

Appendix C. Proof of Lemma 7.1 

The goal of this appendix is to prove Lemma 7.1. The commutation formula (7.1) has 
already been proved at the beginning of section 7. Thus, it only remains to prove the 
commutation formula (7.2). Recalling (3.25), 

□0 = -dQ{do(p) + A(p + n-^Vn ■ V0, 

Thus, we have: 

[□,A]0 = [-dodo + n-^Vn-V + A,A](I) (C.l) 
= -[9o9o,A]0+[n-Vn- V,A]0. 

We thus have to calculate the commutators [Oq, A](f) and [n~^Vn • V, A]0. For any tensor 
U tangent to S^, we denote by VqU the projection of DqU to E^. We have the following 
commutator formula for any vectorfield U tangent to E^: 

[Vb, VoPa = hcVcUa - n-^VbUVoUa + {n'^kabV - n-^kbcVan + Roabc)Uc, (C.2) 

while for a scalar 0, the commutator formula reduces to: 

[Vfe, Vo]0 = kbcVc<l) - n" Vf,n<9o0. (C.3) 

Using the commutator formulas (C.2) and (C.3) and the fact that [do, A](f) — [Vo, V"] Va0+ 
V"[Vo, Va]0, we obtain: 

[do, A](/) = -2A;"''V„V60 + 2n-^VbnVb{do(l)) + n-^Andocf) - 2n- V„nA;"''V6(/), (C.4) 

where we used the constraint equation (2.2) and the fact that, in view of the Einstein 
equations and the symmetries of R, we have: 

g"''Roa6c — 0. 
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Differentiating tlie commutator formula (C.4) with respect to do and using the commutator 
formulas (C.2) and (C.3), we obtain: 

= -2r'VaV5(9o0) + 2n-'VbnVb{do{do<p)) + (-2Vor' + VaV,0 
+ (2n-Vfe((9on) - lOr^'n" Van)Vb((9o0) + (n^^An + 2n-'^\Vn\'^)do{do(i)) 
+ (2r"Roacfe + ^k^'Vckab - 2n- VanVoA;"'' + 2k''''n-^Va{don) + Ak'"'ka,n-^V aU 
^-2\k\^n-^Vbn - 2k''^n-'^V andon)Vb<t> 
+{n-^A{don) - 4r^n-^VaV6n + 2n-^VbnVb{don))do(l). 

Together with the commutator formula (C.4) applied to do(f), we obtain: 

[dodo,A]<f> (C.5) 
= [9o,A]9o0 + 5o([ao,A]0) 

= -Ak'^'VaVbidocf)) + 4n-'VbnVb{do{do(f))) + (-2Vor'' + Ak'^'^k, ") V^V^^ 
+{2n-^Vb{don) - 12r^n-V„n)V6(9o0) + (2n-^An + 2n-^\Vn\^)do(do(l)) 
+{2k'"'Roacb + 2k'"'Vekab ' 2/1" V^nVoA;"" + 2k'"'n-^V a{don) + Ak'^'^k^bU-^V 
+2\k\^n~^Vbn - 2k"-^n~'^W andQn)Wb(p 
+(n-^A(don) - Ak^'^n-^VaVbU + 271'"^^ bnVb{don))do(t). 

We also compute the commutator [n"^ VnV, A]0: 

[n- VnV, A]0 

= -A{n-^Vbn)Vb(t> - ^ ain'^^ bn)V b(t> + n-^Vbn[Vb, A]0 
= -n-^Vb{An)Vb(t> - rT^ [A, VbJnVb^ + V„nV„ V^nVfe^ 
+n-^V6nVanV„V6(/) - V„V6nV„V60 + VfenfVfe, A](/). 

Now, we have the following commutator formula: 

[Vfe, A](/) = i^ft =Vc0 = (Rfeoo ' + kbak'^)Vc(t>, (C.6) 

where we used the Gauss equation for R, the Einstein equations for R and the maximal 
foliation assumption. Thus, we obtain: 



[n- VnV, A](j) (C.7) 
= {-n-^WaS/bn + n-'^VbnVan)VaVb(t> + (-n" ^^(An) + n-'^VauVaVbU 

+2(R600a + haka '')n-^^ aUjV b(t> ■ 
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Finally, (C.l), (C.5) and (C.7) yield: 

[5o9o,A]0 

= [9o,A]9o0 + 9o([9o,A]0) 

= -4r^V„V,(9o0) +4n-V5nV5(9o(9o0)) 

+ (-2VoA:"^ + Ak'"'h''- n^VaVfen + n'^V buV an)V aS/ b(l> 

+{2n-^Vb{don) - 12rWan)V6((9o0) + {2n-^An + 2n"^| Vnn(9o((9o0) 

+ (2r'^Roacfe + ^k^^Vckab - 2n- VanVoA:"^ + 2r^n- Va(9on) + 4A;""A;c6n- V„n 

+2|A;|^n"V6n - 2k''^n''^V audon - n~^Vb{An) + n^'^V bn 

+2(Rwoa + kbaka'')n~^Van)Vb<i) 

+ (n-^A((9on) - 4r WaV^n + 2n-'^VbnVb{dQn))do(t), 
from which (7.2) easily follows. This concludes the proof of Lemma 7.1. 
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